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ABSTRACT 

There  is  a  need  for  this  country  to  be  able  to  assure  itself  that 
an  international  ban  on  nuclear  testing  can  be  made  effective.  This 
thesis  studies  one  small  portion  of  that  problem  -  to  wit,  the  procedure 
for  classifying  an  underwater  detonation  which  has  vented  and  which  has 
been  detected  by  a  hydroacoustic  network. 

A  MEASURE  OF  EFFECTIVENESS  is  developed,  and  after  a  study  of  some 
parameters  such  as  height  of  search,  search  spacing,  search  width  and 
current  movements  optimizing  search  plans  are  deduced. 

One  of  the  major  difficulties  in  such  a  quantitative  analysis  is 
the  non-deterministic  nature  of  many  of  the  parameters.  The  ideal 
search  spacing  was  found  to  depend  directly  on  the  surface  radioactive 
pool  diameter.  In  order  to  develop  plans  this  essential  factor  must  be 
known.  Unfortunately  this  depends  on  many  variables  over  which  the 
searcher  has  no  control  or  knowledge.  An  averaging  procedure  is  adopted 
and  a  factor  established  to  compare  spacings  and  determine  which  meets 

the  criterion  dictated  by  the  MEASURE  OF  EFFECTIVENESS. 
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1 .  Introduction. 

In  October  1960  a  technical  symposium  was  held  by  the  Advanced 
Research  Projects  Agency  (ARPA)  acting  in  behalf  of  the  Department  of 
Defense  for  the  purpose  of  informing  those  in  attendance  of  the  scope, 
status  and  planning  underway  for  Project  VELA.  Project  VELA  consists  of 
research  and  development  that  is  aimed  at  the  improvement  of  our  capa- 
bility of  detecting  and  identifying  nuclear  detonations  that  are  con- 
tained below  the  surface  of  the  earth,  or  underwater  or  that  occur  in 
space  above  the  atmosphere. 

It  is  well  recognized  that  it  is  possible  for  a  clever  and  deter- 
mined nation  to  conduct  subsurface  and  high  altitude  nuclear  weapons 
tests  without  a  significant  risk  of  being  detected.  Detection  of 
nuclear  detonations  on  the  surface  or  within  the  atmosphere  is  consid- 
ered well  in  hand  and  not  included  in  VELA.  Ever  since  the  Conference 
of  Experts  held  in  1958  in  Geneva,  international  groups  have  been 
attempting  to  establish  methods  and  network*  which  would  assure  the 
participating  nations  that  agreements  to  ban  nuclear  testing  could  be 
effectively  enforced  and  thereby  be  the  first  step  in  solving  one  of  the 
most  important  problems  of  our  time  -  ARMS  CONTROL. 

The  hope  of  Project  VELA  is  to  improve  our  detection  and  classifi- 
cation capability  by  at  least  one  order  of  magnitude,  and  perhaps  by  as 
much  as  three,  over  a  three  year  span.  Such  a  quantum  jump  is  considered 
feasible  if  sufficient  independent  scientific  interest  can  be  aroused 
while  assigning  a  high  priority  to  the  investigation  in  as  many  of  the 
government  laboratories  which  have  some  measure  to  contribute. 

The  Project  was  subdivided  into  three  major  categories: 
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VELA  HOTEL  -  Detection  of  very  high  altitude  nuclear  detonations  by- 
means  of  detectors  mounted  in  satellites; 

VELA  SIERRA  -  Development  of  techniques  and  equipment  to  detect  and 
identify  high  altitude  nuclear  detonations  by  ground  based  instrumenta- 
tion} 

VELA  UNIFORM  -  Improve,  markedly,  equipment  and  techniques  for  de- 
tecting and  identifying  subsurface  nuclear  explosions  from  a  distance  as 
well  as  to  develop  rapid,  economical  techniques  for  conducting  on-site 
inspections  to  determine  whether  a  suspicious  event  was  actually  a  nuclear 
explosion  or  not. 

Project  VELA  UNIFORM  includes  a  study  of  underwater  detonations  as 
part  of  its  "subsurface"  objective.  The  Naval  Electronics  Laboratory 
(NEL)  at  San  Diego  was  assigned  the  task  of  investigating  the  detection 
of  underwater  explosions  and  the  Naval  Radiological  Defense  Laboratory 
(NRDL)  was  contracted  to  devise  a  searching  procedure  for  obtaining  sam- 
ples to  classify  explosions  detected  by  a  network  devised  by  NEL. 

In  a  CONFIDENTIAL  report  /6/  NEL  outlined  a  distribution  of  hydro- 
acoustic  stations  which  would  be  capable  of  detecting  underwater  blasts 
in  most  areas  of  the  world.  However,  this  network  has  practically  no 
ability  to  classify  explosions  as  nuclear  or  non-nuclear. 

NRDL  has  pursued  the  problem  of  classifying  a  reported  detonation 
by  locating  the  scene  by  means  of  an  appropriate  search  based  on  the  fix 
information  received  from  hydroacoustic  and/or  other  correlated  equip- 
ment, obtaining  samples  if  possible  and  analyzing  pertinent  character- 
istics of  the  area. 

The  classification  problem  splits  into  two  distinct  categories  re- 
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fleeting  two  completely  different  physical  characteristics  of  the  scene. 
In  the  first  case,  the  underwater  nuclear  explosion  vents  to  the  surface 
and  scatters  the  attendant  radioactive  debris  on  the  surface  and  probably 
into  the  air.  In  the  second  oase  the  radiation  is  contained  below  the 
surface  due  to  a  combination  of  depth  of  explosion  and  yield  of  the 
device.  This  thesis  is  devoted  to  analyzing  the  first  of  these  two 
categories • 

2.  Statement  of  the  problem. 

The  objective  of  this  study  is  to  determine  the  optimum  searching 
procedure  for  locating  the  site  of  an  underwater  explosion  which  has 
vented,  in  order  to  classify  it  as  nuclear  or  non-nuclear. 

3.  Means  of  identification. 

The  most  positive  means  of  confirming  that  an  explosion  is  nuclear 
is  by  a  radiochemical  analysis  of  a  sample  of  the  affected  area. 

After  an  explosion  has  been  detected  and  a  radiating  sample  obtained, 
it  remains  to  distinguish  the  origin  of  the  radiation.  In  spite  of  a 
radioactive  sample  a  nuclear  detonation  could  be  disclaimed  and  the 
activity  attributed  to  one  of  the  following: 

a.  Residue  from  previous  explosions, 

b.  Waste  dumping, 

6.  Reactor  discharge, 

d.  Accident. 

Radiochemistry  can  detect  fission  fragments  and  normally  determine 
the  time  of  detonation  by  extrapolating  back  from  ratios  of  certain  iso- 
topes present  in  the  sample.  This  information  correlated  with  seismic, 
hydroacoustic,  air  acoustic,  and/or  electromagnetic  information  of  time 
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and  location  of  detonation  can  readily  refute  a  claim  of  residue  from 
previous  explosions  or  waste  dumping  as  the  source. 

Normal  reactor  discharge  does  not  contain  fission  fragments  and 
would  easily  be  repudiated  as  a  cause  for  contamination.  However,  leak- 
age from  a  reactor  primary  system  could  be  a  source  of  deceiving  activity 
in  the  sample,  but  would  not  explain  the  detected  explosion  at  the  loca- 
tion. 

An  explosion  of  a  reactor  in  water  would  have  attendant  debris,  but 
with  the  burden  of  proof  to  counter  such  a  claim  upon  the  investigator, 
failure  to  find  debris  is  unlikely  to  be  satisfactorily  conclusive  that 
such  an  explosion  did  not  take  place. 

Any  surface  sample  is  subject  to  prior  masking  by  the  dumping  of 
fission  fragments  to  confuse  analysis.  A  determined  effort  to  make 
radiochemical  analysis  of  surface  samples  inconclusive  is  almost  bound 
to  be  successful  in  a  water  contaminated  area  if  the  deceiver  is  willing 
to  invest  enough  effort.  Contaminated  air  samples  from  venting  of  under- 
water bubbles  are  less  vulnerable. 

Underwater  pulsing  of  the  bubble  on  its  route  to  the  surface  leaves 
a  trail  of  discrete  thin  laminae  of  radioactive  material.  To  reduce  the 
possibility  of  confusion  of  the  surface  sample  by  a  masking  attempt  a 
further  effort  should  be  made  to  sample  some  of  these  lower  strata. 

Sampling  procedures  should  recognize  biological  phenomena  which  can 
corroborate  other  evidence.  Certain  algae  have  a  known  affinity  for 
Iodine  131  which  has  a  short  half  life  and  can  be  identified  as  a  fission 
fragment.  Plankton  can  also  produce  evidence  of  an  event.  If  there  are 
land  and  shallow  water  in  the  vicinity,  molluscs  and  tropical  plant  sam- 
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pies  could  produce  valuable  information. 

Of  course,  specifically,  well  equipped  search  vehicles  oould  arrive 
at  reliable  conclusions  at  the  scene  with  pulse  height  analysis  equip- 
ment, air  drop  water  probes  and  well  schooled  personnel.  However,  it  is 
not  realistic  to  imagine  planes  particularly  outfitted  for  this  searching 
chore.  Our  forces  must  be  more  flexible  for  the  large  variety  of  mis- 
sions they  may  be  called  on  to  perform.  It  is  envisioned  that  patrol 
craft  could  be  rigged  with  the  detector  described  in  USNRDl>TR-525/l/  or 
that  bases  from  which  these  planes  operate  could  have  these  meters  avail- 
able. However  the  more  sophisticated  equipment  is  not  likely  to  be  as 
well  distributed. 

Actual  meter  readings  would  be  of  some  assistance  in  correlating 
with  the  time  of  the  explosion  by  determining  decay  rates j  however,  con- 
tinued dilution  due  to  wind  and  sea  movements  will  make  analysis  of  this 
data  much  less  reliable  than  water  samples. 
4..  Detection  capability. 

It  is  worthwhile  to  consider  some  of  the  major  factors  involved  in 
the  transmission  of  sound  waves  in  water. 

The  sound  reaching  a  receiver  experiences  three  kinds  of  degradation: 

a.  Divergence  -  For  a  point  source  in  a  homogeneous  medium  the 
energy  is  dispersed  in  a  sphere  as  it  is  radiated  out  so  that  a  receiver 
perceives  only  an  amount  of  energy  proportional  to  the  percent  of  area 
of  the  expanding  sphere  that  it  can  observe. 

b.  Attenuation  -  There  is  a  damping  of  the  energy  because  the  pres- 
sure is  not  in  phase  with  the  condensation  due  to  viscosity,  heat  conduc- 
tion and  molecular  absorption  (due  to  magnesium  sulphate). 

c.  Anomalies  -  There  are  many  anomalies  which  contribute  to  the 
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peculiarities  of  underwater  sound  transmission  but  the  principal  ones  are 
reflection  and  refraction. 


The  velocity  of  sound  in  a  fluid  (C)  =  J(B) 

1(P 


(P) 

where  B  =  Bulk  Modulus 
P  =  Density 
As  Temperature  goes  up  P  decreases; 
As  Salinity  goes  up  B  increases; 

As  Pressure  increases  P  increases  but  B  increases  more. 
Therefore  the  velocity  increases  as  temperature,  salinity  and  pressure 
increase. 

Pressure  increases  with  depth  but  temperature  normally  decreases,  and 
temperature  reductions  have  an  overriding  effect  on  the  velocity.  In  the 
Atlantic  a  temperature  of  2°C  is  reached  at  about  4-000  feet  and  remains 
fairly  constant  from  then  on.  Salinity  remains  relatively  constant  and 
velocity  of  sound  Increases  with  depth  due  to  the  pressure  increase.  In 
the  Pacific  the  same  conditions  obtain  at  about  2000  feet.  Near  the 
equator  this  condition  of  constancy  is  reached  at  a  lesser  depth. 

Near  the  surface  radical  temperature  changes  are  not  at  all  infre- 
quent. Surface  disturbances  of  storms  mix  the  water  for  a  considerable 
depth  and  sometimes  send  warm  water  down  while  the  colder  water  is 
brought  to  the  surface.  This  produces  temperature  inversions  which  have 
a  drastic  effect  on  sound  transmission.  This  will  be  discussed  in  more 
detail  below. 

As  the  velocity  increases  the  ray  (perpendicular  to  the  wave  front) 
curves  downward.  Conversely,  if  the  velocity  decreases  the  ray  curves 
upward. 
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Consider  the-  following  ray  diagrams  where  the  left  depicts  the  tem- 
perature changes  with  depth  and  the  right  shows  the  ray  path  (Oo  is 
original  direction  of  propagation). 
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Fig.  1.  Effect  of  temperature  changes  with  depth  on  sound  ray  path. 

Sound  channels  occur  when  there  is  a  minimum  ▼elocity  at  some  depth 
with  higher  velocities  above  and  below.  The  region  of  minimum  velocity 
is  called  the  axis.  Sound  originated  on  the  axis  and  propagated  along 
the  axis  will  not  diverge  spherically  but  only  cylindrically  which  means 
that  the  losses  are  cut  down.  Additionally,  there  is  no  major  anomaly 
to  degrade  the  wave. 

a.  Normal  Sound  Channel 
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b.     Pseudo  Sound  Channel 


TEMPERATURE 


RAY  PATH     ^Reflection  from  surface 


D 
E 
P 
T 
H 


Refraction  due  to  temperature  inversion 


c.  At  depth  of  A000  feet  in  Atlantic  and  2000  feet  in  Pacific: 
TEMPERATURE  RAY  PATH 
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Fig.  2.  Sound  channels. 

Sound  sources  do  not  have  to  be  in  the  sound  channel  in  order  for 
sound  to  be  propagated  along  the  channel.  Rays  entering  sound  channels 
at  small  enough  angles  (usually  less  than  15°)  with  the  horizontal  may 
be  entrapped  and  then  act  as  if  they  originated  on  the  axis. 

These  few  facts  may  help  the  reader  to  appreciate  the  following 

comment  paraphrased  from  the  preliminary  NEL  report  /l/x 

Historically  it  is  apparent  that  even  small  explosions  under  ideal 
conditions  can  be  detected  at  ranges  of  many  thousands  of  miles  , 
whereas,  the  most  powerful  explosions  under  other  conditions  may 

1   MILES  will  be  NAUTICAL  MILES  (6087  feet)  throughout  this  thesis. 
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go  undetected.   For  example;  one  pound  of  TNT  exploded  on  the  sound 
axis  in  the  open  ocean  would  be  expected  to  produce  a  stronger  sig- 
nal than  a  100  Megaton  explosion  under  other  extreme  conditions. 
Obviously,  then,  the  hydroacoustic  detection  of  a  nuclear  burst  is 
a  complex  problem  and  a  generalized  statement  of  the  capability  of 
such  a  system  is  not  possible. 

The  interfaces  between  the  air  and  the  water  or  the  land  and  the 
water  provide  a  poor  match  for  transmission  of  sound  waves  (high  coeffi- 
cient of  reflection).  As  a  result  surfaoe  bursts  or  underground  bursts 
are  not  well  detected  through  hydroacoustic  means  just  as  underwater 
bursts  are  not  well  detected  by  seismic  devices.  For  this  same  reason 
stations  are  not  able  to  monitor  inland  seas  and  locations  in  lagoons  or 
with  seamounts  between  the  event  and  the  receiver. 

The  velocity  of  sound  in  water  is  about  4-800  feet  a  second  and  it 
is  possible  that  some  stations  will  receive  a  signal  as  long  as  one  hour 
after  detonation  (less  than  3000  miles).  It  will  take  at  least  two  more 
hours  to  correlate  station  information,  evaluate,  and  call  an  alert. 

There  are  many  sources  of  confusion  in  detecting.  Hydroacoustic 
stations  will  be  unable  to  distinguish  between  chemical  explosions  and 
atomic  explosions.  Earthquakes  and  underwater  volcanic  disruptions  when 
detected  will  be  a  concern  and  frequently  non-discernible  from  man  made 
explosions.  The  splash  of  missiles  in  deep  water,  ship  propulsion  noises 
and  environmental  sounds  characteristic  of  the  ocean  can  be  distinguished. 

In  the  Joint  Committee  Hearings  /8/  Dr.  H.  Brown  of  Lawrence  Radia- 
tion Laboratory  estimated  that  there  were  5000  continental  earthquakes 
per  year  that  could  be  suspected  of  being  nuclear  explosions  of  between 
one  and  five  kilotons  and  that  only  a  few  percent  could  be  identified  as 
earthquakes  by  present  criteria.  About  100  disturbances  would  be  sus- 
pected of  being  nuclear  explosions  of  five  kilotons  or  more.  Since  most 
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of  the  continental  seismic  false  alarms  are  caused  by  these  earthquakes 
it  seems  unlikely  that  the  number  of  hydroacoustic  false  alarms  due  to 
natural  causes  would  be  this  order  of  magnitude.  Earthquake  disturbances 
are  not  likely  to  be  transmitted  with  too  much  energy  through  the  inter- 
face with  the  water.  Volcanic  disruptions  below  the  water  may  be  mis- 
leading, but  these  are  not  near  as  numerous  as  continental  earthquakes. 
The  principal  cause  of  false  alarm  in  the  hydroacoustic  network  is  most 
likely  to  be  chemical  explosions  transmitted  by  the  Sound  Channel  mecha- 
nism, 
5.  Alternative  courses  of  action. 

With  our  objective  clear  in  mind  -  to  determine  whether  a  reported 
-explosion  was  nuclear  or  non-nuclear  -  our  problem  is  reduced  to  search- 
ing the  area  in  question  to  determine  if  there  is  evidence  of  nuclear 
activity.  Only  by  not  obtaining  evidence  of  nuclear  activity  can  we 
assume  that  the  disturbance  was  non-nuclear  (ignoring  the  possibility  of 
non-venting  nuclear  explosions  for  this  study),  and  even  then  the  conclu- 
sion is  not  absolute.  The  certainty  with  which  a  conclusion  can  be 
reached  depends  on  the  thoroughness  of  the  search. 

The  courses  of  action  that  are  open  revolve  about  the  approach  to 
the  development  of  a  search  plan. 

a.  Should  the  plan  be  based  upon  a  random  search  approach,  consid- 
ering the  target  to  be  equally  likely  in  seme  large  area  prescribed  by 
hydroacoustic  fix  information?  This  would  permit  parallel  sweeps  through 
the  area  in  an  orderly  fashion  by  several  aircraft  with  no  danger  of  in- 
terfering with  each  other.  It  would  assure  a  thorough  coverage  of  the 
area  because  navigation  is  simple  and  reliable, 
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b.  Should  the  plan  attempt  to  optimise  search  effort  with  respect 
to  time?  This  would  permit  target  distribution  to  be  any  shape  that  the 
physical  situation  dictates.  The  pattern  would  be  formed  to  search  the 
areas  of  greatest  probability  continuously  with  the  distribution  of  the 
fix  being  the  "a  priori"  function  and  an  "a  posteriori"  function  formed 
after  each  search.  This  approach  would  be  best  for  the  following  con- 
siderations t 

(1  )  Search  craft  may  have  to  be  recalled  prior  to  completion  of 
the  plan  or  there  may  be  a  cost  or  utility  consideration. 

(2)  International  ground  rules  may  limit  time  to  search  a  given 
area  or  it  may  be  desired  to  spend  as  little  time  as  possible  in  a  given 
area  due  to  political  or  other  considerations. 

(3)  There  is  the  ever  present  marginal  case  where  just  enough 
radioactivity  is  present  at  one  time  in  sufficient  intensity  (because  of 
yield,  depth  of  detonation,  and/or  dispersion)  to  be  detected,  but  which 
will  not  be  there  at  some  later  time. 

A  disadvantage  to  this  approach  to  the  search  plan  is  that,  unless 
target  location  distribution  is  in  fact  uniform,  the  navigation  problems 
reduce  the  reliability  of  conforming  to  a  given  pattern,  and  holidays 
are  left  in  the  coverage.  The  probability  of  detection  may  in  fact  be 
less  than  a  random  pattern,  although  the  same  coverage  effort  is  invested, 
because  of  those  holidays.  Secondly,  there  is  a  peculiarity  in  the 
radioactive  pool  contour  under  many  conditions  which  reduces  the  proba- 
bility of  detection  at  certain  times.  The  character  of  many  pools  is  as 
follows  t 

Upon  initial  venting  a  high  intensity  area  is  established  for  a 
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certain  distance  around  "ground  zero."  After  same  tine  passes  many  of 
the  fission  fragments  sink,  and  fringe  area  activity  is  diluted  beyond 
recognition.  This  decreases  the  contour  of  a  given  low  isodose  line. 
Then  the  particles  that  are  suspended  near  the  center  spread  out  with 
water  and  wind  movements  and  the  low  isodose  contours  spread  out  again, 
even  beyond  the  original  area. 

It  is  impossible  to  fill  in  details  of  "how  much?"  or  "how  long?" 
or  "what  intensity?"  because  all  answers  depend  on  the  physical  condi- 
tions of  venting  and  water  movement.  The  point  of  mentioning  it  here  is 
to  acknowledge  that  it  could  affect  the  capability  of  our  searching 
effort. 
6.  Measure  of  effectiveness. 

This  situation  is  different  from  many  military  operations  where  the 
accomplishment  of  the  objective  is  imperative,  regardless  of  the  cost. 
The  search  would  be  conducted  during  peacetime  and  there  may  be  concern 
for  the  expense  involved,  particularly  if  many  false  alarms  are  encountered. 

There  is  always  a  possibility  of  search  vehicles  having  to  be  recalled 
because  of  more  urgent  commitments,  logistic  considerations  or  forbidding 
weather.  As  previously  mentioned,  international  agreements  may  limit 
searching  time. 

Finally,  there  is  a  matter  of  using  consistent  logic  in  approaching 
the  problem.  As  presented  so  far,  and  as  the  physical  situation  dictates, 
the  problem  admits  the  probablistic  nature  of  the  fix  information.  Again 
ignoring  the  non-eventing  situation,  there  can  never  be  an  absolute  assur- 
ance that  the  explosion  was  non-nuclear,  but  the  probability  of  this  fact 
approaches  unity  as  the  search  effort  approaches  an  extremely  large  value. 
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Therefore,  the  search  will  terminate,  If  a  pool  Is  not  located,  at  some 
time  before  the  evidence  is  absolutely  conclusive  that  no  pool  exists. 
This  will  always  be  the  case  and  the  amount  of  search  effort,  or  synony- 
mously -  time,  is  inherent  in  the  problem  and  cannot  be  ignored • 

under  these  considerations  the  MEASURE  OF  EFFECTIVENESS  is  to  maxi- 
mize the  probability  of  detection  per  time  of  search  effort. 
7.  List  of  variables  to  be  considered. 

a.  Height  of  search.  This  is  discussed  in  detail  in  Appendix  I* 
The  conclusion  is  to  fly  as  low  as  operationally  feasible.  A  plot  of 
intensity  versus  height  in  Appendix  I  shows  the  loss  of  intensity  for 
altitudes  from  0  to  4-0,000  feet. 

b.  Detector  sensitivity.  Appendix  III  is  based  on  a  detector  which 
is  capable  of  distinguishing  0.2  microroentgens  at  the  meter.  In  order 
to  meet  this  requirement  the  meter  envisioned  will  be  collimatedto  reduce 
the  effect  of  cosmic  radiation. 

c.  Search  spacing.  In  Appendix  II  the  search  theory  involved  is 
discussed,  and  it  is  generally  concluded  that  the  spacing  depends  upon 
the  search  width  (which  is  the  effective  visibility  in  the  definite  range 
law)  and  the  standard  deviation  of  the  location  errors. 

d.  Search  width.  The  radioactive  pool  on  the  surface  of  the  water 
coupled  with  the  width  of  the  detector  view  determine  the  search  width. 
This  is  covered  in  detail  in  Appendix  III. 

e.  Target  motion.  The  movement  of  the  radioactive  pool,  which  is 
the  object  of  the  search,  will  depend  principally  on  currents.  Appendix 
IV  investigates  the  variability  of  the  water  movement  and  concludes  that 
the  speed  of  movement  is  predictable  but  the  direction  is  random. 

13 


f .  Fix  distribution. 

g.  Probability  of  detection.  In  Appendix  III  the  probability  of 
detection  is  computed  for  the  ideal  spacing  -  given  the  search  width. 
Since  the  search  width  is  not  known  in  a  realistic  situation  an  aver- 
aging procedure  must  be  resorted  to.  This  is  developed  in  paragraph  11 
below. 

h.  Search  aircraft  characteristics. 

i.  Navigation  errors  in  locating  the  center  of  the  fix, 

j.  Time  late  at  the  scene. 

8.  Hrdroacoustic  fix  distribution. 

The  location  of  an  explosion  will  be  estimated  by  correlation  be- 
tween hydroacoustic  detection  stations.  Their  estimates  will  depend  on 
physical  condtions  prevailing,  but  for  this  study  the  errors  have  been 
classed  as  having  a  standard  deviation  of  2.5  miles  or  25  miles.  The 
purpose  of  assuming  deviations  an  order  of  magnitude  apart  is  to  see  if 
the  search  plan  is  sensitive  to  such  extremes.  The  distribution  is  as- 
sumed normal  and  the  probability  is  represented  byr 

f (x,y)  =  f  (r)  =  J_  EXP(-  r2) 

27T<r;-   (  W1) 

On  the  basis  of  this  assumption  Appendix  II  develops  the  ideal 
search  spacing  as  0.75  (EC")*  where  E  is  the  Effective  Visibility  and  <T 
is  the  standard  deviation. 

9.  Navigation  error. 

Where  there  are  no  visual  navigational  aids  available  an  aircraft 
can  expect  to  arrive  at  a  given  location  in  accordance  with  some  proba- 
bility distribution.  According  to  Bowditch  /u/  celestial  errors  in  air 
navigation  of  five  to  ten  miles  are  considered  normal  for  favorable  con- 
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ditions  (pg  676  paragraph  2806).  Electronic  navigational  aids,  where 
available,  yield  more  accurate  positioning.  In  this  model  it  is  assumed 
that  the  error  due  to  navigation  is  normally  distributed  with  a  standard 
deviation  of  five  miles.  The  variance  is  additive  to  the  hydroacoustic 
fix  variance  since  they  are  both  normal  and  independent.  When  the  fix 
standard  deviation  is  2.5  miles  the  total  location  standard  deviation  is 
5.6  miles.  When  the  fix  deviation  is  25  miles  the  location  deviation  is 
25.5  miles. 

These  figures  demonstrate  that  the  navigation  error  has  a  substan- 
tial effect  on  the  problem  when  the  fix  error  is  small,  but  it  is  rather 
inconsequential  when  the  fix  error  is  large. 
10.  Search  patterns. 

Appendix  IV  discusses  currents  in  detail  and  for  this  study  the  set 
is  considered  equally  likely  in  any  direction  at  a  given  location  and 
time.  The  drift,  on  the  other  hand,  can  be  predicted  reasonably  well. 
Koopman  /2/  in  sections  1 .6  and  7.3  discusses  a  non-uniform  distribution 
of  targets.  When  the  fix  deviation  is  2.5  miles  the  situation  falls 
closely  into  the  category  of  a  circular  distribution  of  the  target  with 
movement  at  a  known  speed  in  a  random  direction.  The  theory  is  to  search 
on  the  peak  of  the  distribution,  then  outside,  then  inside,  ...  etc. 
This  pattern  will  give  the  maximum  probability  of  detection  per  unit  time 
for  locating  a  moving  target  whose  initial  fix  is  not  definite.  See 
Appendix  V  for  details  of  this  plan. 

When  the  location  deviation  is  of  the  order  of  25  miles  or  more  and 
the  time  late  is  less  than  two  days  the  target  can  be  treated  as  station- 
ary and  the  search  pattern  that  applies  is  described  by  Koopmen  /2/  in 
paragraph  7.3.2.  A  set  of  expanding  squares  -  squares  of  uniform  cover- 
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age  -  are  superimposed  in  the  search  pattern.  See  Appendix  VI  for  de- 
tails of  this  plan, 
11 •  Search  width. 

This  one  parameter  remains  to  be  determined.  Appendix  III  discusses 
the  geometry  for  determining  the  search  width  and  it  concludes  that  it  is 
the  width  of  a  uniform  pool  of  minimum  intensity  minus  the  detector  view. 
Unfortunately,  there  is  no  way  to  determine  the  pool  size.  There  are 
really  not  enough  statistical  data  available  to  generally  characterise 
pool  contours  for  various  depths,  yields,  sea  and  wind  conditions  and  as 
a  function  of  time. 

a.  Large  location  error. 

As  a  matter  of  interest  the  following  tables  have  been  computed 
to  consider  the  effect  of  using  a  search  width  for  one  size  pool  when 
the  pool  is  actually  some  other  size.  These  calculations  are  based  on 
the  estimated  probability  of  detection  given  by  Koopman  /2/. 

P  =  1  -  (1  -taw )  exp  -nw 

(  d       r 

n  =  Number  of  searches  made  (one  for  this  table). 
w  =  Search  width, 
s  =  Spacing  of  pattern. 

As  concluded  in  Appendix  II  the  spacing  is  .75(20*  and  for  this 
problem  W  =  E.  The  standard  deviation  considered  is  25»5  miles. 
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TABIE  1.     Analysis  of  probability  of  detection  for  given  search  spacings, 

Search  spacing   (S)  =  A.T1  mil—  -  Computed  for  a  2  mile  diameter  pool 

Pool  Diameter  Search  Width        W  Probability  of 

(Pool  of  uniform  minimum  intensity)  (W  =  E)             S  Detection 

2  1.3               .302  0.036 

4.  3.3              .766  .18 

6  5.3  1.23  .348 

8  7.3  1.692  .505 

10  9.3  2.16  .634 

20  19.3  4.48  .939 

30  29.3  6.8  .991 

Search  spacing   (S)  =  6.86  -  Computed  for  a  4  mile  diameter  pool 

2  1.3               .1896  .015 

4  3.3               .481  .085 

6  5.3              .7725  .18 

8  7.3  1.062  .287 

10  9.3  1.352  .394 

20  19.3  2.813  .7714 

30  29.3  4.27  .9262 

Search  spacing   (S)  =  8.7  -  Computed  for  a  6  mile  diameter  pool 

2  '1.3               .1495  .01 

4  3.3              .3795  .055 

6  5.3              .609  .125 

8  7.3               .839  .205 

10  9.3  1.068  .289 

20  19.3  2.22  .652 

30  29.3  3.375  .8513 
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Search  Width 

W 

s 

Probability  of 
Detection 

1.3 

.1273 

.0075 

3.3 

.324 

.04 

5.3 

.52 

.095 

7.3 

.715 

.161 

9.3 

.912 

.23 

19.3 

1.393 

.564 

29.3 

2.875 

.78 

Search  spacing   (S)  =  10.2  -  Computed  for  an  8  mile  diameter  pool 
Pool  Diameter 

2 

4 
6 
6 

10 
20 
30 
Search  spacing   (S)  =  11.5  -  Computed  for  a  10  mile  diameter  pool 

2  1.3  .113  .008 

4 

6 

8 

10 

20 

30 

Search  spacing   (S)  =  16.6  -  Computed 

2 

4 
6 

8 
10 
20 
30 
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3.3 

.287 

.034 

5.3 

.-461 

.076 

7.3 

.634 

.132 

9.3 

.808 

.191 

19.3 

1.678 

.5 

29.3 

2.55 

.723 

a  20  mile  diameter 

pool 

1.3 

.0783 

.005 

3.3 

.1987 

.0175 

5.3 

.319 

.04 

7.3 

.44 

.073 

9.3 

.56 

.107 

19.3 

1.163 

.325 

29.3 

1.764 

.525 

Search  spacing  (S)  =  20.4  -  Computed  for  a  30  mile  diameter  pool 


Pool  Diameter 

• 

Search  Width 

W 

s 

Probability  of 
Detection 

2 

1.3 

.0637 

.002 

4 

3.3 

.1618 

.01 

6 

5.3 

.26 

.029 

8 

7.3 

.358 

.05 

10 

9.3 

.456 

.078 

20 

19.3 

.946 

.243 

30 

29.3 

1.436 

.421 

Figure  3  is  a  plot  of  the  probability  of  detection  for  each  of  the 
search  spacings  under  considerations  against  the  widths  of  the  pools.  A 
plot  against  search  widths  would  be  the  same  shape  displaced  to  the  left 
0.7  miles.  Figure  4  utilizes  the  same  data  but  uses  SPACING  vice  POOL 
DIAMETER  as  the  abscissa. 

In  order  to  consider  this  nebulous  problem  more  quantitatively, 
assume  that  any  search  width  from  0  to  30  miles  is  equally  likely  and 
investigate  the  mean  probability  of  detection  for  each  search  spacing 
over  the  span  of  search  widths. 


P=   /  p(w)f (w)dw  =    /jl-d+W)  exp(-W)]dw. 
•'O  y»   L    S)     SJJ30 


(l-eacp(-W)  -  H  exp(-W) 
L      S)   S     S\ 


W+S  exp(-W)  +  S  exp(-W)  (W+1  ) 
S)        S)  (S  J 


30 


P  =  J.  jWJ 
30  U 

P  n  _l]  30  +  S  «rcp(-20)  +  S  exp(-20)  (£0+1 )  -  2sl 
30L  S)        S)  (  S      J 

The  Measure  of  Effectiveness  is  again  called  upon  to  decide  which 
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spacing  affords  the  best  approach.  It  is  desired  to  select  the  spacing 
which  yields  the  maximum  probability  per  unit  time  of  search  effort.  The 
spacing  (S)  is  inversely  proportional  to  the  time  of  search  effort. 
Hence,  by  maximizing  pS  the  criterion  previously  established  will  be 
satisfied. 

TABLE  2.  Comparison  of  spacings  using  the  MEASURE  OF  EFFECTIVENESS  (LARGE 
LOCATION  ERROR. 


Associated 
Bool  Diameter 

Spacing 
(S) 

2 

4.31 

4 

6.86 

6 

8.7 

8 

10.2 

10 

11.5 

20 

16.6 

30 

20.4 

Aye. 

Irob.  of  Detect. 

Measure  of 
Effectiveness    (T5S ) 

.71355 

3.07 

.561 H5 

3.855 

.47058 

£.10 

.409 

4.17 

.3633 

4.175 

.239 

3.965 

.183 

3.73 

From  the  above  table  it  can  be  seen  that  pS  is  maximized  around  a 
spacing  of  11.5  miles. 

By  differentiating  the  expression  for  p*S  with  respect  to  S  and 

solving  for  the  maximum  the  following  is  obtained: 

d(g5)=  dTs+S2  exp(-20)4s£  exp(-20)(2P+1 )-  SJH 
dS   dS  L  30     S  )  30     S  )  (S     1 5  J 

d(T3S)=0=l+exp(-?0)+  1  S  exp(-20)+20  exp(-2P.)+2  exp(-20)+_S  exp(-20) 
dS  .        S  )  15       S  )  S      S  )       S  )  15     S 

-2S 
15 

The  following  transcendental  results: 
S  In  (2S2+45S+450)  =  30 


(2S1U5S 
(     2S21-, 


22 


S  =  11  satisfies  this  expression. 

Substituting  S  =  11  into  the  expression  for(jJS) 
pS  =  4.18  (5  s  0.38) 

As  a  result  of  the  above  -  based  on  the  very  stringent  assumption 
that  search  widths  will  vary  uniformly  from  0  to  30  miles  (implying  pools 
from  0  to  30.7  miles  in  diameter)  -  it  is  concluded  that  a  search  spacing 
of  11  miles  maximizes  the  ratio  of  probability  of  detection  to  time  of 
search  effort.  This  corresponds  to  anticipating  an  average  search  width 
of  8.44  miles  and  a  pool  diameter  of  9.14  miles  of  uniform  minimum 
intensity. 

If  the  assumption  had  been  that  pool  widths  varied  uniformly  from  0 
to  50  miles  in  diameter  the  conclusion  would  have  been  that  a  spacing  of 
25  miles  is  ideal.  It  is  apparent  that  the  conclusion  is  very  sensitive 
to  this  assumption.  Statistical  data  are  required  to  provide  a  firmer 
base  for  this  conclusion,  but  until  data  becomes  available  the  assumption 
as  stated  is  considered  reasonable. 

b.  Small  Location  Error. 

When  the  standard  deviation  of  the  location  error  is  only  5.6 
miles  the  effect  of  the  current  is  relatively  significant.  After  a  delay 
of  24  hours  the  peak  of  the  distribution  would  have  moved  out  to  a  circle 
of  12  miles  radius  (assuming  a  drift  of  0.5  knots).  The  search  path  is 
along  the  peak  at  the  time  of  arrival  (modifying  spiral  to  an  expanding 
square).  After  the  first  path  the, "a  posteriori"  probability  becomes 
bimodal  with  a  circular  peak  inside  and  one  outside.  The  pattern  con- 
tinues to  search  these  maxima. 

It  is  beyond  the  scope  of  this  thesis  to  develop  a  probability  of 
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detection  statement  for  this  situation,  but  it  is  worth  considering  the 
tendency  of  the  wa  posteriori"  distribution.  As  each  peak  is  searched 
the  surface  of  probability  becomes  flatter  until  after  one  pass  of  the 
area  the  distribution  approaches  a  uniform  one.  As  the  first  search  is 
being  made,  however,  the  situation  is  not  much  different  from  that  de- 
veloped in  the  case  of  the  larger  deviation.  In  the  latter  case  the 
probability  peak  was  continuously  being  pushed  outward.  In  the  situation 
of  the  moving  target  the  circular  peak  is  being  split  into  two  concentric 
rings  enclosing  the  original  peak.  Each  resulting  ring  is  searched,  and 
the  peak  is  pushed  both  inward  and  outward  as  the  search  progresses. 

In  order  to  get  some  feel  for  the  problem  (which  is  more  difficult 
mathematically  than  physically)  the  same  expression  for  the  probability 
of  detection  is  used  as  a  foundation  for  investigating  the  search  width. 

From  a  geometry  standpoint,  after  an  area  of  a  given  peripheral 
size  is  searched,  0.75  would  receive  the  same  intensity  of  search  in  the 
current  problem  as  was  received  in  the  stationary,  large  deviation  prob- 
lem. The  remaining  quarter  would  be  only  half  covered  but  with  twice 
the  effort. 
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Fig.  5  Comparison  of  search  patterns. 
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Under  these  very  rough  approximations  and  the  assumptions  of  the 
previous  discussion  a  table  similar  to  Table  2  can  be  computed. 


Associated 
Pool  Diameter 

Spacing 

(s) 

2 

2.02 

4 

3.22 

6 

4.08 

8 

4.79 

10 

5.41 

20 

7.79 

30 

9.61 

Average  Prob.  of 
Detection  (v) 

Measure  of 
Effectiveness  (ES) 

.767 

1.55 

.70 

2.25 

.654 

2.67 

.615 

2.95 

.582 

3.15 

.472 

3.68 

.409 

3.93 

TABIE  3.  Comparison  of  spacings  using  the  MEASURE  OF  EFFECTIVENESS  (SMALL 
LOCATION  ERROR). 

This  indicates  that  the  largest  spacing  would  maximize  the  approxi- 
mate probability  of  detection  per  time  of  search  effort  provided  the 
actual  probability  of  detection  is  at  least  proportional  to  the  above 
estimation  and  that  pool  sizes  are  equally  likely  in  size  from  0  to  30 
miles. 

After  one  pass  of  the  area  using  the  doubly  expanding  square  and  a 
spacing  of  9.61  miles  the  target  location  distribution  approaches  a  UNI- 
FORM one.  It  is  obviously  not  UNIFORM,  but  the  peaks  have  been  flattened 
and  multiplied  and  the  continued  movement  of  the  current  throughout  the 
search  builds  up  the  probability  of  the  target  being  in  a  path  already 
swept.  The  end  result  is  that  the  unit  mass  is  best  considered  to  be 
evenly  spread  over  the  area. 

As  a  result  of  this  reasoning  a  parallel  sweep  search  will  maximize 
the  probability  of  detection  per  time  of  search  effort  after  the  first 
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pass.  The  expression  for  probability  of  detection  then  becomes: 

P  =  W  f or  S>Wj  P  =  1  for  S£V. 
S 

< 

Following  the  same  approach  as  in  subparagraph  11 .a.  it  is  desired 
to  determine  the  spacing  that  meets  the  established  measure  of  effectire- 
ness  -  again  under  the  assumption  of  equally  likely  search  widths  (pool 
diameters)  over  the  range  of  0  to  30  miles. 

po 

P  =  J.  /   w  d*   for  s>w 
30/   S 
y0 


=  _L 
30 


r  o3° 
2S 


=  15 

s 


We  are  interested  in  pS.  Sinoe  p  cannot  exceed  1  there  is  no  sense 
in  making  S  smaller  than  1 5  because  p  would  remain  constant  as  S  was  made 
smaller  and  pS  would  therefore  decrease.  Obviously  the  maximum  is  reached 
when  S  =  W.  Therefore  this  search  should  be  conducted  at  a  spacing  of  15 
miles. 

For  the  second  pass  the  search  area  should  be  a  square  whose  dimen- 
sions are  dictated  by: 

a.  Three  times  the  location  standard  deviation  to  provide  an  almost 
certainty  that  the  pool  was  included  in  this  circle  at  the  time  of  the 
fix,  plus, 

'  b.  the  product  of  the  drift  and  the  time  that  would  elapse  between 
the  time  of  the  fix  and  the  time  of  completion  of  sweep  of  the  area  on 
this  second  pass. 

The  sweeps  should  start  along  the  side  perpendicular  to  the  most 
frequent  direction  of  the  current  at  15  mile  spacing. 
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12.  Aircraft  capabilities. 

It  is  assumed  that  an  aircraft  will  be  available  for  the  search 
equipped  with  the  proper  survey  meter  and  water  sample  floats  capable  of 
the  following: 

a.  Maximum  range  to  scene  1000  miles  with  minimum  of  3  hours  on 
station  available. 

b.  Search  speed  180  knots. 

13.  Search  procedure. 

This  problem  deals  with  an  underwater  burst  which  vents  to  the  at- 
mosphere. All  of  the  effort  of  the  study  has  been  devoted  to  an  analysis 
of  activity  at  the  scene,  but  the  study  has  presupposed  that  an  air 
sampling  barrier  would  be  established  downwind  from  the  blast.  Any 
samples  obtained  from  this  endeavor  would  not  be  subject  to  masking  and 
quite  informative  about  the  circumstances  of  the  origin  of  the  blast. 

As  previously  mentioned,  correlation  of  hydroacoustic  information 
may  take  several  hours.  It  is  reasonable  to  assume  that  as  much  as  2A 
hours  may  pass  before  search  craft  could  arrive  at  the  scene.  The  appro- 
priate search  plan  should  be  followed  depending  on  the  physical  situation. 
If  a  pool  or  suspected  pool  is  located  an  attempt  should  be  made  to  locate 
the  area  of  highest  intensity  (with  due  regard  for  safety  of  the  crew). 
Intensity  readings  should  be  taken  and  correlated  with  position  and  time 
as  accurate  as  possible.  Relative  positions  within  the  pool  are  more 
important  than  exact  geographical  locations.  At  the  point  of  highest 
intensity  water  sampling  floats  should  be  dropped  for  later  recovery  by 
surface  craft  or  helicopters. 

The  search  aircraft  should  remain  in  the  area  as  long  as  safety  and 
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operational  requirements  permit  so  it  can  direct  other  searching  activity 
for  recovering  water  sample  floats. 

Frequent  intensity  readings  may  provide  sufficient  information  to 
establish  decay  curves  to  extrapolate  back  to  the  time  of  detonation. 
Great  accuracy  would  be  necessary  in  correlating  the  movement  of  the  pool, 
height  of  aircraft,  location  of  aircraft  over  some  relative  point  of  the 
pool,  etc.  to  permit  much  of  an  estimate  of  the  decay  rate.  The  best 
procedure  would  be  to  keep  recording  the  intensity  of  the  area  of  highest 
intensity  to  provide  some  semblance  of  a  reference.  This  would  corrobo- 
rate hydroacoustic  information. 

Normal  search  procedures  should  be  followed  in  marking  datum  and 
search  area,  and  making  reports. 
14.  Conclusions. 

The  problem  was  to  determine  the  optimal  searching  procedure  for 
locating  the  site  of  an  underwater  explosion  which  had  vented,  in  order 
to  classify  it  as  nuclear  or  non-nuclear. 

A  hydroacoustic  network  is  anticipated  to  be  the  detecting  agent. 
Such  a  network  has  no  ability  to  distinguish  an  underwater  nuclear  explo- 
sion from  a  non-nuclear  one.  Classification  must  be  accomplished  by  a 
thorough  investigation  of  the  scene  for  evidence  of  nuclear  activity.  If 
no  such  evidence  is  found  a  probability  statement  can  be  made  of  the 
nature  of  the  explosion,  and  upon  this  information  a  decision  to  terminate 
or  continue  the  search  can  be1 reached. 

A  MEASURE  OF  EFFECTIVENESS  was  established  to  be: 

Maximize  the  probability  of  detection  per  time  of  search  effort. 

Several  parameters  were  investigated  to  appreciate  their  effect  on 
the  problem. 
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a.  Height. 

It  was  determined  that  the  maximum  intensity  reading  at  the 
meter  occurs  on  the  surface.  As  altitude  is  increased  the  intensity 
decreases  up  to  an  altitude  above  20,000  feet.  After  that  the  increase 
is  gradual  and  almost  negligible.  Of  course  at  such  a  height  the  in- 
terest is  only  academic  since  very  few  waves  have  sufficient  energy  to 
survive  that  thickness  of  air. 

b.  Search  Spacing. 

The  development  for  search  spacing  as  presented  by  Koopman  /2/ 
was  expanded  and  adopted  for  this  study, 

c.  Search  Width. 

The  effective  search  width  was  determined  to  depend  on  the  radio- 
active pool  diameter  which  was  an  unknown  quantity.  Pool  diameters  depend 
on  yield,  depth  of  detonation,  and  physical  characteristics  of  the  water 
and  the  exploded  device.  In  order  to  discuss  the  problem  quantitatively 
an  averaging  procedure  was  used,  assuming  any  pool  diameter  to  be  equally 
likely  between  0  and  30  miles.  This  is  a  very  brash  assumption,  and  the 
results  of  the  study  depend  heavily  on  this  weighting  function.  Unfortu- 
nately, little  statistical  data  are  available  which  might  indicate  the 
distribution  of  this  function.  This  assumed  UNIFORM  distribution  of 
weights  provides  am  unexpected  conclusion  when  the  criterion  is  the  MEASURE 
OF  EFFECTIVENESS  mentioned  above.  If  the  weight  function  favored  some 
particular  pool  diameter  it  would  be  reasonable  to  adopt  the  most  effect- 
ive spacing  for  that  diameter;  however  if  there  is  no  bias  toward  a  par- 
ticular width  the  conclusion  is  not  so  evident.  Various  spacings  were 
concluded  to  be  optimum  depending  on  the  standard  deviation  of  the  dis- 
tributions considered. 
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d.  Location  Error. 

There  are  two  categories  of  location  error  considered.  Large 
errors  with  standard  deviations  above  25  miles  dictated  a  search  plan 
that  ignored  current  drift.  For  smaller  deviations  the  drift  had  to  be 
reckoned  with  in  designing  the  search  plan.  The  location  error  consid- 
ered the  effect  of  both  the  initial  hydroacoustic  fix  error  and  naviga- 
tional errors  of  aircraft  arriving  at  a  given  spot  in  the  open  ocean. 

When  current  movement  was  considered  (with  small  standard  deviations 
of  location)  the  direction  was  assumed  to  be  random,  but  the  drift  can  be 
predicted  within  an  acceptable  margin. 

e.  Search  Plans. 

When  the  standard  deviation  of  the  location  distribution  is 
large  the  squares  of  uniform  coverage  maximize  the  probability  of  detec- 
tion per  time  of  search  effort.  On  the  other  hand  if  the  deviation  is 
small, a  retiring  search  maximizes  this  criterion  for  the  first  pass.  For 
subsequent  passes  of  the  small  deviation  problem  the  "a  posteriori" 
rippled  distribution  approaches  a  UNIFORM^and  a  parallel  sweep  search  is 
considered  best. 

Finally,  although  this  comment  is  beyond  the  scope  of  the  stated 
problem,  it  seems  worth  mentioning  in  passing  that  very  little  special 
equipment  is  required  to  perform  this  function  of  classifying  detected 
disturbances.  Water  sample  floats  which  can  be  dropped  from  an  aircraft, 
sample  the  water  and  seal  themselves  -  to  be  recovered  by  surface  craft 
or  helicopters,  and  the  survey  meters  described  by  NRDL  /l/  are  the  only 
non-standard  items  needed.  These  items  can  be  based  ashore  and  on 
carriers  and  need  not  be  on  aircraft  until  a  plane  is  assigned  a  particu- 
lar task  of  classifying. 
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APPENDIX  I 
HEIGHT  OF  SURVEYING  AIRCRAFT 

In  order  to  determine  an  appropriate  search  plan  the  search  height 
must  be  determined.  According  to  reference  /l/  the  meter  will  be  con- 
structed with  a  shield  forming  a  cone  about  the  axis  to  prevent  excessive 
interference  from  cosmic  radiation.  This  shielding  forms  a  central  angle 
of  150°. 

As  the  search  craft  increases  its  height  the  area  seen  by  the  meter 
will  increase  and  one  would  expeot  this  to  increase  the  intensity  at  the 
instrument.  However,  the  attenuation  due  to  the  atmosphere  will  decrease 
the  intensity  at  the  meter. 

Model  assumptions t 

1 .  The  pool  is  UNIFORM  and  INFINITE  in  diameter. 

2.  Gamma  energy  at  the  surface  is  1  MEV  and  the  "build-up" 
factor  is  not  considered, 

3»  Absorption  is  linear  varying  at  a  constant  rate  with  dis- 
tance from  the  source. 

4..  Height  of  meter  is  much  greater  than  the  diameter  of  the 
meter  opening, 

I  =  Intensity  of  radiation  at  the  meter, 
A  a  Area  of  meter  opening. 

I  =  Intensity  of  point  sources  at  the  surface, 
h  =  Height  of  instrument. 

u  a  linear  coefficient  of  absorption  in  air  (2.265  x  10  /FT). 
x  =  Distance  of  meter  from  radiating  source. 

r  =  Distance  of  source  along  water  from  point  on  surface  directly 
below  aircraft. 
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0  =  Angle  of  rotation  of  wrw  on  surface, 

0  =  Angle  that  point  radiating  source  makes  with  meter. 

w  =  Solid  angle  subtended  by  meter  opening  at  radiating  source. 

The  portion  of  the  radiation  that  arrires  at  the  meter  from  each 

point  source  is  equal  to  the  ratio  of  the  solid  angle  (w)  subtended 

divided  by  U   . 

W_   A  cos  0 
U    ~  x2  U 

This  assumes  that  h  is  much  greater  than  the  diameter  of  A;  there- 
fore the  projection  of  the  area  of  the  meter  opening  on  a  plane  tangent 
to  a  sphere  expanded  about  the  point  source  is  approximately  equal  to 
the  projection  onto  the  sphere, 

{project**^ 


Ail,    \ 
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Fig,  7  Parameters  of  integration, 
h  tan  0     /-27T 


I  = 


I     A  cos  0    exp[-ux7  r  dr  dO 
^ffJZ 
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cos  0  -  h. 

X 


r  =  (x2-h2£ 
dr  = 


(x2-h^>£ 


(1)  1  = 


h  sec  0 


L  Ah  expf-uxjdx 
2  "  air 


It  would  be  of  interest  to  find  a  height  which  would  maximize  the 
intensity  at  the  meter  if  such  exists. 

The  integrand  in  (1  )  is  continuous  except  at  o  and  the  derivative 
can  be  investigated  by  applying  Leibnitz's  Rule: 


(2)  dj  =  I^A 


dh 


h  exp(-u  h  sec  0)  sec  0  -  h  exp(^u  h) 
b/  sec2  0  h^ 


h  sec  0 

eap(-wx)dx 


"j 


Consider  the  last  integral  in  (2), 

h  sec  0 

exp.(-ux)dx= 


/*        ~]h  sec  0 
-exp(-*ux)-u  j   exj)(-ux)dxj 


["-|h  sec  0 
-exp(-ux  )-u  (in  x-ux+ujx^-ail+a^  . . 
x  2*2!  3*31  4*4! 

—  h 

This  term  by  term  integration  is  justified  by  virtue  of  uniform 

convergence  of  the  infinite  series  expansion. 

Substituting  this  in  (2)'. 

(3)  dl  =  -In  sec  0  +  u  h(seo  0-1  )  -  u2h2  (sec2  0-1  ) 
dh  U 


+  Ji2h2(sec3  0-1  )  -  u^-bA  (sec^-1  ). 
18  96 

Sec  0  =  3.8637 
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dl  was  programmed  on  the  1604  computer  for  a  large  number  of  values 
dh 

of  h.  The  results  shoved  that  the  derivative  was  always  negative,  be- 
coming smaller  in  absolute  value  as  height  increased.  This  indicates 
that  a  plot  of  the  intensity  would  be  an  exponential  decay  with  its  maxi- 
mum value  at  h  =  0. 

Assumption  3  in  the  previous  model  is  not  realistic  because  the 
attenuation  varies  with  density.  Since  air  density  varies  with  height 
the  coefficient  of  absorption  must  also  vary  with  height.  The  mass  co- 
efficient of  absorption  is  equal  to  the  linear  coefficient  divided  by 
the  density  and  it  does  not  vary  with  the  chemical  or  physical  condition 
of  the  medium.  For  air  this  quantity  is  0.0274  anr/gm. 

According  to  the  Barometer  equation  derived  in  STATISTICAL 

MECHANICS i 

P  =  P0  exp(-3g£) 
KT 

P  =  Density  at  height  of  instrument. 

m  =  Mass  of  a  molecule  of  the  medium. 

g  =  Acceleration  due  to  gravity. 

h  =  Height  of  the  instrument. 

K  =  Boltzman's  constant. 

T  =  Temperature  (degrees  KELVIN). 

PQ=  Atmospheric  density  at  sea  level. 

u»  =  Mass  coefficient  of  absorption. 

x  =  Slant  range  from  radiating  source  to  meter. 

0  -   Limit  of  shielding  (75°). 

Returning  to  equation  (1 )  and  substituting! 
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«xp  [j-tDcJ   =  aocpij-u'RcJ   =  exp  UufxP0  exp(-mgh  )  J 

u  KT     J 


I  = 


vh  860  0 


I»  Ah  exp  Ut^xP,.,  expC-ggh)  dx 
2         x^L  °  KT  ;J 


(4.)    1  dl  =  expJ^£»P0  exp(-mgh)h  sec  0J  -  ffiE[^a,P0  exp(-mgh)hj 
B  dh       hsec0 


KT 


KT 


h  see  0  — 

h  mg  exp(-mgh)ufPff  ex^f-u'PQ  eocp(-mgh)x  dx 


KT  KT 


KT 


(5) 


s-h  sec  0 
+  exp  [Ha1  P^>  exp(-mgh)xjdx 

The  last  integral  of  equation   (4.)  becomes i 

h  sec 


B  =  IoA 
2 


h  sec  0 
-  exp(-ax)! 

'h  h 


/ 


exp(-ex)  dx 

x  Where  a=ufP0  exp(-mgh) 

KT  y 


(6) 


The  first  part  of   (5)  balances  the  first  two  terms  of   CO. 

/h  sec  0 
exp(-ax)dx 
x 
n 

In  sec  0  -  ah  (sec  0-1  )  +  a2h2  (sec2  0-1  )... 

2*2!  -» 


=  a(h  mg  -1  ) 
KT 


Consider  the  following* 

"a"  is  positive  for  all  values  of  h; 
The  infinite  series  is  always  positive; 

h  mg  -  1  is  0  when  h  =  KT 

KT  mg 

Therefore  there  is  a  maximum  or  minimum  or  flex  at  h  =  KT 

mg 

d2I  =  B J  2n  sec  0-a  h(sec  0-1  )  +  a2h2   (sec2  0-1  )>  >  . 

dhZ    h  =  KT         X  2«2! 


n« 
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Again  "a*,  B  and  the  series  are  positive.  The  conclusion  is  that 

h  =  KT  is  a  minimum  point, 
mg 

Equation  (6)  was  programmed  out  to  30  terms  for  values  of  height 
from  -100  to  40,000  feet  and  several  sensitivity  checks  were  made. 

The  temperature  is  recognised  as  a  variable  over  the  range  investi- 
gated and  it  was  varied  from  300°K  (81  °P)  to  273°K  (32°F).  The  values 
varied  slightly  but  the  structure  of  the  curve  remained  the  same.  A  few 
points  were  checked  with  the  temperature  at  243°K  (-4,0°F)  and  again  the 
height  of  the  minimum  decreased  slightly  bui  the  structure  of  the  curve 
remained  the  same.  This  covered  the  range  that  average  temperatures 
could  reach  and  the  minimum  varied  from  21100  feet  to  26070  feet» 

Variations  in  surface  air  density  (P0)  were  looked  at  but  their 
effect  were  even  less  significant. 

Variations  in  "g"  were  not  looked  at,  but  they  are  obviously  of  no 
concern. 

Lastly,  the  evaluation  was  made  for  10  terms  of  the  series  and  then 
for  30  terms  of  the  series.  In  both  cases  the  values  were  the  same  out 
to  eight  significant  figures. 

As  a  final  investigation  the  expression  for  the  intensity  was  pro- 
grammed and  curves  obtained  are  presented  in  figure  8.  It  is  satisfying 
that  the  same  MINIMUM  points  were  reached  by  two  independent  approaches 
and  two  different  programs. 

The  development  was  as  follows: 

h  sec  0 
I  =   /         IQA  h  exp(-ax)dx 
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-  InAh  -exp(-ax) 
2     L.X 


h  sec  0 


-  ~  I/v^fa)  «cp(-h  sec  0  a)-exp(-ha)  +a  In  sec  0-a  hfsec  0-1  )... 
2[h  sec  0  h  J     L  --J 

=  lo&\  h~l _Vhaf3a  sec  0+O-l\ha(sec  0-1 ) 

2    [\     sec  0)      '  2!'; 

V-  1    ^h2a2(sec2  0-1  )...-/+  1      -    1     V  an(sech0-1  ).../ 
V  2-2!  "J  (-h^hl'+Th+TT!']  -1 

I  =  2aA|.74.1+ha(-la  sec  0+d-l  )  ha  (sec  0-1  ) ) J 

In  fig.  8  the  bracketed  quantity  is  plotted  as  f (h)  -  the  "Intensity 
factor". 
CONCLUSIONS 

1 .  The  greatest  intensity  at  the  meter  will  be  experienced  when  the 
height  is  zero. 

2.  The  minimum  point  of  intensity  is  little  more  than  academic 
interest  since  the  assumption  of  a  UNIFORM  source  of  diameter  at  least 
equal  to  the  field  of  the  meter  becomes  a  stringent  condition  at  25,000 
feet.  The  pool  would  hare  to  be  30  miles  in  diameter.  More  directly, 
this  problem  is  interested  in  maximum  intensities. 

3.  Operational  conditions  will  dictate  the  height  of  search  with 
the  aim  of  flying  as  low  as  feasible  on  a  prolonged  mission.  Figure  8 
demonstrates  the  degradation  of  intensity  at  the  meter  (and  therefore 
the  decrease  in  search  width)  for  the  height  of  search  selected.  It  is 
not  possible  to  predict  the  decrease  in  probability  of  detection  which 
will  be  experienced  as  height  increases  because  the  physical  conditions 
of  the  pool  of  radiation  are  not  known.  Since  the  meter  has  a  threshold 
intensity  which  must  be  reached  before  a  signal  can  be  recognized,  flying 
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at  a  greater  altitude  than  operationally  necessary  may  degenerate  the 
intensity  at  the  meter  so  that  the  threshold  is  not  reached  -  whereas 
it  may  have  been  reached  at  the  lower  level. 

A»     At  the  altitudes  that  are  operationally  feasible  the  curve  will 
give  a  reliable  indication  of  Intensity  factor  (f (h))  for  any  atmospheric 
conditions  (neglecting  the  build-up  factor  effects )  since  it  is  insensi- 
tive to  temperature  changes  below  3000  feet. 
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APPENDIX  II 
SEARCH  THEORY 
This  appendix  is  based  on  Koopman's  Search  and  Screening  (OEG  Report 
No.  56)  /2/.  An  attempt  is  made  to  present  his  development  in  more  de- 
tail as  it  applies  to  this  thesis', 

A.  Development  of  probability  of  detection  of  a  target  of  lateral 
range  X  from  a  straight  aircraft  track t 
(1  )      p(X)  =  1  -  expr-0.092/E^2]] 

1.  Consider  the  case  of  "continuous  looking"  where  "s"  is  the 
instantaneous  probability  density  of  detecting  a  target.  If  s  remains 
constant  throughout  a  search  p(t)  =  l-exp(-st)  where  p(t)  is  the  prob- 
ability of  detection  during  time  t. 

If  s  is  a  changing  density  because  of  relative  movement  with  respect 
to  the  observer,  and  is  time  -  dependent,  the  expression  is  more  compli- 
cated. The  instantaneous  density  s  is  inherently  dependent  on  range  and 
the  fact  that  it  varies  with  time  is  conveyed  by  the  symbol  s^. 

p(t)  =  1    -  expj-   Jst  dt] 


db&f*r 


Fig.  9  Target's  relative  track. 


The  target  is  at  (x,y)  at  time  t  and  x  =  x(t),  y  =  y(t).  Let  the 
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initial  position  at  t  =  t»  be  XQ  =  X(tf),  70  =  y(t»),  and  the  final  posi- 
tion at  t  =  t"  is  ^  =  X(t"),  y,  =  y(t")  with  the  target  tracing  a  track 

along  C  in  figure  9  above. 

s  =s[jX2(t)4r2(t))i]::St 

The  probability  of  detection  along  path  C  is: 


£  r  *  s([x2(t)4y2(ti]^)dtj 


p.  =  1  -  exp. 


The  integral  is  a  line  integral  along  C;  if  "w"  is  the  relative 

speed  and  nl"  is  the  arc  length  along  C  from  XQ,  Y0,  then 

pc  =  1  -  expL-  js(Qc2(t)+y2(t)]^)  §1  dt  =  dl 

c  w  w 

Following  Koopman's  notation  the  exponent  is  called  the  sighting 
potential  F[Cj.  This  quantity  possesses  the  property  of  additivity  so 
that  where  C  =  0,  +  C2,  then  F[c]  =  F^CJ)  +  f£c^| 

In  the  case  where  X  is  the  lateral  range,  X  is  constant  and  y  =  wt 
(w  is  a  constant), 

fJcJ  =      r  *'  .{  kWt2j*>4t  =1      f  T"  ■(  kVJ^My 

v  r 

According  to  the  inverse  cube  law  of  sighting  where  h  <<  f 
(2)  s  =  k_h. 


'tG]  "  ¥  f  '"(i^p/a 


yf 


=  Lh 

w 


y1 


"  *$\P  <**+/$ "  x^s5^5^! 
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Ob<?<2.<*^3  r 


y"«x-x 


Fig.  10  Detection  at  fixed  speed  and  course. 

When  the  sweep  is  "relatively"  infinite  ©f  =  ©"  =  90°  and  - 

When  the  observer  and  the  target  are  on  straight  courses  at  constant 
speeds  for  a  long  time  before  and  after  their  closest  point  of  approach, 
the  probability  of  detection  is  a  function  of  the  lateral  range  X. 

^<*K  /-eyf  [-^-*2 

2.  Parallel  sweeps  - 

A  target  is  presumed  to  be  at  rest  on  the  ocean  in  an  unknown 
position,  all  equal  areas  having  the  same  chance  of  containing  it.  A 
search  is  made  along  a  large  ("infinite")  number  of  parallel  lines  a 


common  distance  apart.  This  distanoe  is  referred  to  as  Sveep  Spacing  (s) 
and  the  problem  is  to  determine  the  probability  of  detection  of  a  target, 
P(s),  by  this  configuration.  Other  questions  concerning  probability  of 
detection,  if  target  speed  and  course  are  known  or  have  a  distribution 
of  their  own,  can  follow  in  development. 

Consider  the  lateral  ranges  of  the  target  from  various  sweeps 


(*.tf> 


X 


l*z    ■    - 


-IS  -5      o  s 

Fig.  11   Parallel  sweeps. 

For  paths  to  the  left  the  distances  are  X,  X  +  s,  X  +  2s...  and  for 
those  to  the  right  s  -X,  2s  -X,  3s  -X...  Both  cases  can  be  combined 
in  the  absolute  expression 

lateral  range  =  X  -  ns     where  O  <t   X  <  S 

n  =  0,  ±1,  ±2,  ±3,  ... 

The  probability  of  detection  by  the  n**1  lateral  sweep  when  the  tar- 
get is  at  (X,  y)  ist 

The  probability  of  no  detection  by  the  n^h  sweep  is  1  -  p^;  and  no 
detection  by  any  sweep  is  the  infinite  product  "TT  (1  -  Pq  )  for  all  values 
of  n.  The  probability  that  at  least  one  sweep  detects  a  target  at 
(X,  y)  is  ^ 

co     ^x,*j  =  /- jr^rL-^iv-i*!^/-  exf  L-0(\s2 


<X3 


(5)  where      (p^  «\ 


.    £     fCU-nSl);       o<y<s 
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Consider  the  lateral  distance  of  the  target  X  as  having  a  uniform 
distribution  between  0  and  s.  P(s)  Is  the  average  of  P(X,  s)  over  all 
values  of  X  in  the  interval.  The  expected  or  average  value  of  the 
probability  of  detection  for  a  given  s  can  be  expressed  as  a  function 
of  s. 

(6)  PCS) *£%C\s>£ji**  iS^O-^H'  l-$(^)^jdK 

This  is  the  solution  to  the  problem  but  it  remains  to  express  it 
more  quantitatively. 

The  effective  visibility  (E)  is  defined  as  half  that  sweep  spacing 
for  which  the  probability  of  detection  by  parallel  sweeps  is  \\   i.e. 
P(2E)  =  \. 

In  the  definite  range  law  detection  occurs  if  and  only  if,  the  tar- 
get happens  to  be  within  the  definite  range  R  of  either  of  two  adjacent 

sweeps.  The  chance  for  this  is  2jl  =  W  when  s  >2R  =  W,  and  unity  when 

s    s 

s  ^.  W.  Prom  this  it  follows  that  E  =  W  (-since  W  =  .5 J  JL  =  *5j  W  =  E). 

s      2E 

For  the  inverse  cube  law  when  h  <^<  r  we  get  from  equation  (3 ) » 

F(X)  =  2  h  kt  and  from  equation  (5)r  0(1,   s)  =  2  h  k 
VTX*  T 

This  can  be  restated  as  fS   (X,  I )  =  2  k  h  J  csc2^  because  of  the 

V     s2  s 

following  developments 

a.  The  infinite  product  expression  for  sin  S  is 
sin  8  =  ZTT    (I-.    J—  ,  ) 

b.  Take  In  of  each  side  - 


In  sin  I  =  -&  £  t     2*?(\-  ftjj^) 

U6 


c.  Differentiate  both  sides  - 


sing     ^        ^.-iU!    ^rrw   v.     h*.^/ 


d.  Which  may  be  rewritten: 

cot  2  =  4-  -  S 


Note:  This  is  the  Fourier  expression  for  cot  2  and  could  have  been  used 
as  the  starting  point  of  this  development. 
e.  Break  up  d.  and  solve  for  A  and  B  by  partial  fractions: 


-  JL  ( '  JL-  t  -£-  Y 

From  partial  fractions  = 

A(nff+  »)  +  B(nlT-  2)  -  22 

A2  -  B2  =  22 

A  =  2  +  B 

AnTT+  BnfT=  0 

2  +  B  +  B  =  0  .        B  =  -1 

•    » 


(7)     Co+  €  ~~  ^   '    fc^   [ 


A  =  2  +  B  =  1 


(s)  cot  £  -  j  -+  5      J-_ 

(9)      COt^rs     £        -J—^ 

The  only  way  to  see  the  above  steps  is  to  substitute  values  for  n  in 

(7)  and  see  that  you  get  two  terms  for  each  n  which  can  be  represented  by 

the  expanded  limits  of  the  expression  in  (8).  Note  also  that  when  n  =  0 

in  (8)  the  missing  term  is  1  which  has  thoughtfully  been  provided. 

2 

f.  Differentiate  (9): 
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do)     cscz^  =    S^  (^-HirV 


g.     Let  3  =  T[Xand  rearrange   (10) 

8  cxi? 


TT 


C^C 


*    *      £ 


-4»-5> 


(V~^ 


u 


(11)  Then     (p(X^)  »    £ 


iA 


c^ 


OcT-        J&3  (|  %-hO 


)J<* 


Inserting  this  in  (6)i 

There  is  one  more  bit  of  wizardry  which  Koopman  applies: 


A 


Let  B  = 


ur  S  ■ 


Tl  > 


0(B) 


Differentiate  with  respect  to  B 
01  (B)  = 


Jo*0 

esc2©  d©  =  -d(cot  ©) 

when  ©  =  0    cot  ©  =  fl>« 

©  =  ^  cot  ©  =  0 
2 

Let  cot2©  =  X2 

d  cot  ©  =  dx 

01(B)=   -^   ltf(?&>C°  9.y.f(- &*?■•)  di 


Integrate  0*  (B)  with  respect  to  (B)  from  0  to  B'  (a  particular  value 


of  B) 


4* 


Let  X2  =  B  2  X  dx  =  dB 

X     =  Ef£  when  B  =  0    X  =  Oj 

,  when  B  =  B1 ,  X  =   (B1  ft 

(12)  0(B)  =  ^  J  *  $*p(-  tf-j  1/ 

(B1  ft  -  X1   =  a  particular  value  of  X. 

By  definition   (12)  is  the  error  function   (erf)  based  on  the  normal 

distribution. 

Hence  P(S)  =  0(B)  =  erf   (B)^  =  erf  ^bkhft 

s  V  \r' 

For  the  inverse  cube  law  the  effective  sweep  width  W  =  2(2  khYz 

\  W"/ 

(13)  P(s)  =  erf  W  (m  i 

s  2 

P(2E)  =  £  =  erf  (ff)fr  V 

Since  erf  0.477  =0.5,  then  (7H^V  =  .477 

and  W  =  1 .076E,  hk  =  .0^6E2 

Substituting  this  in  (13) 
P(s)  =  erf  (0.954E) 

8 

It  is  of  particular  interest  at  this  time  to  note  that  the  proba- 
bility of  detection  is  not  particularly  sensitive  to  the  law  by  which 
the  spacing  (s)  is  determined  -  particularly  at  lower  probabilities. 
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I.o 


f(4) 


R  a  ^  e 

Lav 


n  =  1 
s 

Fig.  12  Probabilities  vith  parallel  sweeps. 

Recall  the  following  equations: 

p(X)  =  1  -  exp[^F(X)J  for  lateral  range  distribution 
(l4)p(X)  =  1  -  expf-2hk  1  for  inverse  cube  law  ; 

LvrFJ  ; 

W  =  1 .076E  from  the  error  function  development  • 
(15)  hk  =  0.04.6E2  from  the  error  function  development 

vr 

Substituting  (15)  into  (U): 
P(X)  =  1  -  expf-0.092/E)27 

Which  is  equation  (1 ),  the  immediate  object  of  the  above  expansions. 
Before  proceeding,  let  us  review  the  principal  assumptions  upon 
which  the  "error  function"  development  is  based. 

a.  Target  at  rest  in  an  unknown  position  with  lateral  distribution 
UNIFORM  between  sweeps. 

b.  Observer  is  on  a  straight  course  at  a  constant  speed  for  a  "long" 
time  before  and  after  the  closest  point  of  approach  to  the  target.  (Sweep 
is  "relatively  infinite"). 

c.  Inverse  cube  law  of  sighting  applies  and  h«r.   ("h"  is  the 
vertical  height  of  the  observer  and  "r"  is  the  horizontal  distance  from 

50 


observer  to  center  of  the  target.) 

d.  Sweeps  are  independent. 
B.  SPACING  (S). 

Consider  the  SQUARE  CF  UNIFORM  COVERAGE  search  pattern  (section 
7.3.2  of  reference  2). 


■*s 


4-S 


Fig.  13     Square  of  Uniform  Coverage. 
Let  the  normal  density  function 

(16)  <f  OS  -h)  -  f  (r)  :  ^  e<*  (-  ^0 

represent  the  probability  that  the  target  at  the  time  of  the  fix  be  in 
the  small  region  dx  dy  at  the  point  (x,y),  a  distance  of  r  from  the  origin 
and  where  <T  is  the  standard  deviation  of  fix  accuracy.  The  first  problem 
is  to  find  a  spacing  so  that  on  the  initial  square  the  probability  of 
detection  per  unit  time  shall  be  a  maximum. 

As  was  indicated  before,  the  spacing  is  not  sensitive  to  the  law  of 
detection,  and  the  inverse  cube  law  development  will  be  followed  so  that 
the  probability  of  detection  of  a  target  of  lateral  range  x  from  a 
straight  aircraft  track  is 


(17)  p(x)  =  1   -  expUO.092 


<KJ 
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Reflect  for  a  moment  on  our  problem  and  the  assumptions  of  the  pre- 
vious section.  The  same  error  funotion  result  is  being  used. 

1 .  The  lateral  distribution  is  certainly  no  longer  uniform  between 
sweeps;  however  when  the  standard  deviation  is  large  compared  to  the 
spaoing  the  approximation  of  uniformity  between  sweeps  is  not  very  in- 
accurate. 

2.  The  sweeps  are  not  "infinite"  in  any  one  given  direction,  but 
except  in  the  center  the  length  of  sweep  is  relatively  large,  and  again 
the  practical  approximation  is  acceptable. 

3.  The  inverse  cube  law  is  assumed  in  the  development  with  h<Cr. 
First  of  all,  as  already  stated,  the  spacing  is  not  sensitive  to  the  law. 
Caution  must  be  used  in  proceeding,  however,  because  although  the  distance 
to  the  center  of  the  detected  configuration  can  normally  be  expected  to 

be  greater  than  the  height  of  the  aircraft,  in  the  case  of  an  extremely 
weak  source  approaching  the  minimum  detection  capability  of  the  meter  - 
when  directly  over  the  pool  -  the  range  to  the  center  would  be  less  than 
the  height  and  the  probability  of  detection  would  be  overestimated.  The 
instantaneous  probability  (s)  would  equal  k  h   a^06  ^h  and  the  error 

would  increase  as  the  ratio  of  height  to  horizontal  distance  increased. 

U»     Sweeps  are  still  considered  independent. 

The  probability  of  the  target  lying  on  the  strip  parallel  to  the  y 
axis  between  x  and  x  +  dx  is  found  by  integrating  equation  (16)  over  all 
values  of  y  and  obtaining  the  marginal: 

The  peak  of  the  distribution  is  located  at  the  center.  If  an 
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indefinite  straight  flight  is  made  along  the  y  axis  and  has  failed  to 
detect  the  target  the  distribution  is  altered  and  the  differential  co- 
efficient of  (18)  no  longer  represents  the  lateral  density  of  the  tar- 
gets. By  Bayes  theorem  the  "a  posteriori"  probability  considers  the  "a 
priori"  probability  of  (18)  above  and  the  "productive"  probability  of 
not  detecting  the  target  on  a  sweep  through  the  center: 

l-p(*)-c  £Xf  H-  °'oiil(^TS  (See  ec*uation  (17)) 

Bayes  Theorem  states: 

Let  the  union  of  events  B.  be  the  sample  space  where  Bj 
(i  =  1 ,  2,  ...,  k)  are  mutually  exclusive  and  exhaustive. 

Let  event  C  also  be  defined  on  the  same  sample  space 

S     C]      '  ii    iii iTJrF-    i- 1—  ■—  -i rmi,     i  _ 

in  I 
In  our  problem  let: 

C  =  not  detecting  target 

Bj=  target  located  in  strip  (x,  x+dx) 

The  denominator  of  Bayes  theorem  can  be  expressed  in  the  continuous 

case  as    e^-a 

J-  o-"       QJ    \       %  '    Q  where  y  represents  any  strip . 

This  integral  is  just  a  number,  say  K,  which  actually  normalizes  the 

numerator  so  that  instead  of  speaking  of  the  probability  of  the  event  Bj 

given  event  C  we  speak  of  the  density  function 

-  so  that  the  new  density  function  is  proportional  to  P(CJB. )f(B. )  and 
this  translates  to 
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(i9)  Li- "*>U<*  «  *■'**[-&- '•"HWJ 

K  * J  — 

In  words  -  the  probability  of  not  detecting  the  target  given  that  it 
is  in  a  strip  (x,  x+dx)  times  the  probability  that  it  is  in  a  strip 
(x,  x+dx)  is  proportional  to  the  new  density  of  the  probability  that  the 
target  is  in  a  strip  (x,  x+dx)  given  that  it  was  not  located  on  the  first 
sweep. 

To  find  the  maximum  of  the  new  distribution  differentiate  (19)  and 
solve  for  the  x  which  makes  the  expression  zero. 

So  the  maximum  probability  target  position  is  no  longer  humped  at 
zero.  Double  peaks  now  occur  with  a  depression  at  the  origin.  However 
the  distribution  is  skewed,  and  the  distance  between  the  first  and  second 
tracks,  D,  should  be  slightly  greater  than  .65 (EG")"?  to  cover  a  larger 
volume  of  probability  on  the  second  sweep. 

The  probability  of  not  detecting  on  the  second  sweep  given  that  the 
target  is  in  a  sirip^x,  x+dx) times  the  probability  that  it  is  in  a  strip 
(x,  x+dx)  is  proportional  to  the  new  density  of  the  probability  that  the 
target  is  in  a  strip  (x,  x+dx)  given  that  it  was  not  located  on  the  second 
sweep.  This  follows  from  the  same  reasoning  used  in  equation  (19) 

Integrating  this  function  for  all  x  and  solving  for  a  D  which  wo\xld 
minimize  the  remaining  function  we  would  be  maximizing  the  probability 
of  detection  on  the  second  pass.  Koopman  presents  the  solution,  arrived 
at  by  numerical  integration,  in  section  7.3.2  of  SEARCH  AND  SCREENING  /2/ 
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as  0.75 (E<f)2.  This  spacing  is  then  used  for  searching  about  a  point  of 

fix. 

C.  Search  Patterns 

1 .  Square  search  for  a  stationary  target  (Square  of  Uniform  Cover- 
age"). 

The  details  of  the  theory  are  described  very  explicitly  in 
SEARCH  AND  SCREENING  /2/.  This  search  assumes  a  circular  distribution 
of  target  position  as  a  result  of  a  probablistic  fix.  This  pattern 
attempts  to  give  a  high  probability  of  detection  per  unit  time  at  the 
beginning  of  the  search  when  covering  the  area  where  the  target  is  most 
likely  to  be. 

2.  Retiring  square  search  for  a  moving  target  (Section  7.3.3  of 
SEARCH  AND  SCREENING  /2/). 

In  this  development  it  is  assumed  that  the  direction  of  movement 
is  not  known  but  uniformly  distributed  (i.e.,  equally  likely  in  any  direc- 
tion), while  the  speed  of  movement  is  reasonably  well  established. 

The  theory  is  to  search  the  peak  of  the  distribution  at  the  time 
of  arrival.  It  has  moved  outward  from  the  fix  a  radius  of  Dt  where  U  is 
the  speed  of  target  movement  and  t  is  time  late.  During  the  search  the 
target  continues  to  move  and  the  pattern  is  to  search  the  peak,  then  out- 
side, then  inside,  then  outside,  ...  etc.  This  will  give  the  maximum 
probability  of  detection  per  unit  time  of  search  effort. 

It  is  recognized  that  both  of  these  patterns  present  operational 
problems  in  implementing.  Expanding  searches  are  always  frowned  upon 
because  of  the  navigational  difficulty  in  making  legs  precise.  This  is 
particularly  true  when  spacing  is  small  -  less  than  5  miles. 
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APFENDIX  III 
SEARCH  WIDTH 

USNHDL  Technical  Report  525  /l/  envisionsan  airborne  detection  in- 
strument capable  of  sensing  radioactive  surface  activity  on  the  sea  when 
the  intensity  at  the  meter  is  0.2  microroentgens/hr.  At  a  search  height 
of  500  feet  the  intensity  of  a  broad  source  at  the  surface  would  have  to 
be  0.74.  microroentgens  per  hour.  The  meter  envisioned  would  be  colli - 
mated  about  a  central  angle  of  150°  to  reduce  cosmic  interference. 
A.   For  determination  of  the  effective  search  width  consider  the  fol- 
lowing mod el : 

1 •  Physical  parameters : 

a.  Search  height  -  500  feet. 

b.  Meter  described  above  (i.e.,  shielded  around  arc  of  150°, 
and  capable  of  detecting  0.2  micror/hr). 

c.  Minimum  surface  intensity  present  -  0.74  micror/hr. 

d.  Meter  system  time  constant  about  2  seconds.   (The  technical 
report  anticipates  0.2  sec  for  the  meter  and  it  is  presumed  that  a  re- 
corder, alarm  and  other  auxiliary  equipment  would  make  system  time  con- 
stant an  order  of  magnitude  greater). 

2.  Assumptions: 

a.  Surface  pool  of  radiation  is  circular  and  of  uniform  inten- 
sity. 

b.  The  radiation  on  the  surface  ends  abruptly  at  a  definite 
boundary. 

c.  Definite  Range  Law  of  detection  obtains  (W  =  E). 

Due  to  the  geometry  the  meter  will  see  a  circular  area  about  0.613 
miles  in  diameter  (See  Fig.  H)« 
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Fig.  14.  Meter  geometry  (side  view). 

If  a  circular  area  of  0.74-  micro-r/hr  intensity  and  0.613  miles 
diameter  were  present  in  the  search  area,  the  effective  visibility  (E) 
would  be  zero  since  the  plane  would  have  to  fly  directly  over  the  center 
in  order  to  detect  it. 

If  the  pool  were  larger,  the  search  width  could  be  considered  as  the 
diameter  of  the  pool  minus  the  meter  visibility  (0.613  miles). 
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Fig.  15  Surface  geometry. 

Another  consideration  in  determining  the  effective  visibility  (E)  is 
the  meter  system  response  time.  In  the  USNRDL  report  /1/  the  following 
standard  formula  is  provided. 
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d)  s  =  s0d  -  «p£-t/tJ) 

where  S  =  meter  indication 
SQ=  source  strength 
t  =  time  over  source 
T  =  system  time  constant 

Neglecting  attenuation  and  build-up,  if  the  meter  were  over  the 
source  for  ten  seconds  the  meter  would  indicate  over  99%  of  the  source 
strength. 

For  an  aircraft  traveling  at  180  knots  the  meter  would  have  to  be 
over  the  source  for  0.5  miles.  This  cuts  down  the  effective  diameter 
another  0.1  miles  for  a  pool  of  two  miles  diameter  and  a  lesser  amount 
for  larger  pools. 
OB  =  (0A2-AB2)£ 
OB  =  0.64.6  mi 
EB  =  1 .292 
OF  =  1  mi 
AC  a  0.5  mi 
AF  -   0.3065  mi 
0A  =  0.6935  mi 
C  =  First  meter  indication 
A  =  Last  meter  indication 
CA  =  Time  meter  registers 

over  pool  (10  sec=£mi) 


New  Effective: 
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Fig.  16  System  time  constant  effect. 
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For  this  problem  a  correction  of  -0.7  miles  is  made  to  pool  diameter 
to  obtain  effective  visibility. 

Since  the  size  of  the  pool  is  an  unknown  in  the  search  several  sizes 
will  be  looked  at  to  determine  sensitivity. 

For  a  search  of  a  stationary  object  using  a  Square  of  Uniform  Cover- 
age pattern  assume  a  fix  error  (2  standard  deviations )  to  be  50  miles  and 
an  error  in  navigation  of  10  miles  in  locating  point  of  fix.  The  total 
deviation  would  be  25.5  miles.  For  simplicity  we  will  adopt  the  follow- 
ing expression  for  the  probability  of  detectiont 
(2)  p  =  1-(l+nW\expRiVn    (See  Koopman  /2/  pg  113  and  1U) 

n  =  #  passes  over  area;  W  =  search  width  =  E;  S  =  spacing. 


(miles ) 
Pool  Dia 

E=W 
Effect. 
Visibility  (mi) 

3A(E  StdDev)t= 
Spacing  (miles ) 

W 
S 

}   Prob. 

of  detection 

r 

I   n=1 

n=2 

n=3 

2 

1.3 

4.31 

I  .302 

t 

.024 

.121 

.23 

4 

3.3 

6.86 

.481 

.085 

.25 

.42 

6 

5.3 

8.7 

.609  1 

-. 

.125 

.344 

.545 

« 

7.3 

10.2 

l 

.715  I 

)   .U 

.416 

.632 

10 

9.3 

11.5 

.809; 

j   .193 

.479 

.698 

20 

19.3 

16.6      I 

1.16  | 

.321 

.673 

.857 

30 

29.3 

20.4 

1 .435  I 

.42 

.752 

.928 

TABLE  4.   uniform  square  coverage  probabilities. 

Since  it  is  unlikely  that  data  are  available  for  such  small  intensities 
(0.74  micro  r/hr),  it  is  desirable  to  consider  intensities  that  are  more 
likely  to  be  recorded  and  scale  effects  from  available  search  data. 

Consider  a  circular  pool  of  higher  intensity,  say  0.05  r/hr.  If  the 
attenuating  distance  of  two  configurations  were  equal,  then  for  each  to 
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induce  the  same  intensity  at  the  meter,  the  product  of  each  area  and  sur- 
face intensity  must  equal  the  same  product  of  the  others. 
(3)     J-!  A,^  1?  " z  Lj   =  Minimum  detectable  surface 

6x?Y  t/0~    ff  (O'lobT)-  ,oSflz  intensity  when  spread 

(\    =  U.XT>~  ild"*0  sl  y*c'  uniformly  over  A1 

I2  =  Uniform  intensity  of  new 
configuration 

^_ _____      A_  =  Area  of  new  configuration 
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Fig.  17  Geometry  of  detector  view  of  high  intensities. 

Area  of  radiation  in  view  of  detector  is  the  sum  of  the  segments t 


60 


(5)     Common  chord  yields t        JT      -.      5i  m  -%? 

K.  S>iv>   Jgi, 

The  effective  search  width  =      2./Y?  -f-J,    —  ("  r]  -  A  y  >) 

=   2th, +1 J 

From  CO* 

From  (5):      J[L     Sfvi   -r^    -,    S\  vl  -&- 

(8)      ^  A-dose?,  )-  £\-c*s  6^) 

'3- 

Equations  (7)  and  (8)  are  transcendental  and  can  be  solved  graph* 
ically. 


1  ■ 

From  equation  (7) 

From  equation  (8) 

miles 
Pool  radius 

radl ana 

rgdl ana 

•l 

°fL 

*1 

H 

ALL 

0 

.082 

'" ■ 

U 

0 

.5 

.005 

.008 

1 

.016 

2 

.033 

3 

.049 

4 

.005 

.081 

.065 

5 

.005 

.080 

.081 

10 

.005 

.075 

.163 

.5 

.01 

.015 

1 

.01 

.032 

2 

.01 

.08 

.01 

.066 

3 

.01 

.077 

.01 

.097 

A 

.01 

.073 

.01 

.130 

5 

.01 

.066 

.01 

.162 

10 

.01 

0 

.01          |        .327 

.5 

1 

2 

.015 

.0745 

3 

.015 

.062 

4 

5 

10 

.015 

0 

i i. 

(ContM 

next  page) 
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Pool  radius 


2 


.5 


.5 


.5 


-*i 


.5 


From  equation  (7) 


•02 

*3L 


^02i_ 


.03 


.035 

U021 


.04 

uQL 


^05. 


«2 


.06 


.079 

.065 


.076 
.047 


.072 


.06 


SXl 


From  equation  (8) 


*2 


TABLE  5.  Data  for  graphical  solution  to  transcendental  equations. 
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Fig.  18  is  a  plot  of  these  curves.  The  points  of  intersection  yield 
the  solution  for  0<| ,   ©2  **or  various  pool  sixes. 

Note  the  limiting  case  in  which  the  pool  is  so  large  that  the  arc 
approaches  the  chord  of  the  detector !s  circle.  In  this  situation  all  of 
the  area  is  described  by  the  segment  of  the  detector's  circle. 

A  =  i  r22  (©2  "  8in  °2) 

©2=  0.082  radians. 

The  solutions  for  ©2  range  from  .0675  radians  (3°51 ' )  for  £  mile 
radius  pool  to  0.082  radians  U°42')  for  an  infinite  pool.  With  these 
figures  and  for  the  pools  under  consideration  the  effective  search  width 
varies  from 

2 (pool  radius  +  .306  mi)  for  large  pools  to 

2(0.5  +  .306)  for  the  half  mile  radius  pool. 

Therefore  the  search  width  is  Increased  by  the  diameter  of  the  de- 
tector "look"  at  the  surface  (approximately)  regardless  of  pool  site. 

Remember  that  a  specific  example  of  a  pool  of  0.05  r/hr  with  the 
detector  at  500  feet  altitude  is  the  case  under  consideration. 

There  is  one  very  unrealistic  assumption  in  the  above  discussion. 

The  intensity  at  the  meter  is  not  just  proportional  to  the  area. 

t 
The  distance  of  the  meter  from  the  source  is  of  principal  concern  when 

considering  attenuation.  Ignoring  the  build-up  factor  the  interrelation- 
ship is  as  follows: 


JA  WK1  ^      ***** 


(9) 


.2  micro  r/hr  =  f^^  e^f-^Y)  «-<ircl& 


u 


(The  effect  of  the  Barometer  equation  is  considered  negligible  -  see 
Appendix  I) 

Equation  (9)  would  be  solved  numerically  for  Aj  or  more  meaningfully 
for  the  12  required  in  the  pool  to  register  at  the  meter. 

A  rough  approximation  for  the  effect  of  the  different  distances  of 
the  source  from  the  meter  can  be  obtained  by  comparing  the  intensity  of  a 
particle  at  the  average  distance  of  each  configuration. 

The  average  distance  to  the  particle  in  the  case  of  the  circle  of 
minimum  intensity  is 

jf  rdr  d&  TTr^ 

(For  circle)  Z^   =  0.223 

In  the  case  of  the  segment  the  average  distance  is: 

X2  =  .3178 

I*t  I-|  =  meter  intensity  for  circle  of  minimum  intensity 

I2  =  meter  intensity  for  segment 

L|0  =  surface  intensity  for  circle 

I20  =  surface  intensity  for  segment 
I  =    ToAlo   &X£  C-^l   for  point  source. 

Foi^   I^-K^f  (-^E^JT?,1    1^0  0^\ l   * 

In  order  for  the  meter  intensity  to  be  the  threshold  intensity  in 
the  case  of  both  patterns,  the  surface  intensity  of  the  segment  would 
have  to  be  10.65  the  intensity  of  the  oircular  configuration  to  compen- 
sate for  the  different  mean  distances. 
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In  the  previous  scaling  discussion  the  characteristics  deduced  would 
pertain  to  a  surface  intensity  of  0.5  r/hr  vice  0.05  r/hr. 
The  following  table  will  assist  in  scaling t 
Altitude  =  500  ft;  £,  =  .223 %   u  =  13.788/mi. 


Il0(r/hr) 

A2(mi2) 

Limiting 

eg 

h2(»i) 

X2(mi) 

^oAio 

J20 

W(mi)  = 
Fbol  Size 
plus 
(h2+r2) 

.5 

4.375X10-7 

2.18 

.3064 

.3175 

10.6 

5 

.6129 

.05 

4.375X1 0"6 

4.7 

.3063 

.3175 

10.6 

.5 

.6128 

.005 

4.375X10-5 

10.14 

.3057 

.3167 

10.4 

.05 

.6122 

.0005 

4.375X1 0-4 

21.94 

.3015 

.314 

9.8 

.005 

.608 

.00005 

4.375X10-3 

47.73 

.2805 

.311 

9.1 

.0005 

.587 

.000005 

4.375X10-2 

107.72 

.1809 

.235 

— — , 

1.4 

i . 

.00005 

.487 

If  search  altitude  is  1000  feet  the  following  table  can  be  used: 
Z|  =  .446 


Il0(r/hr) 

A^Oni2) 

Uniting 
°2 

hgdni) 

X^ini) 

^cAio 

X20 

W(mi)= 
Pool  + 

.05 

1.003X10-3 

10.5 

.6101 

.635 

39 

1.95 

"^U22 

.005 

1. 003x1 0~2 

22.5 

.6001 

.628 

34.6 

.173 

i 

1.21 

.0005 

1. 003x1 0"1 

48.9 

.559 

.599 

19.9 

.00995 

1.17 

.00005 

1.003 

111.2 

.346 

U565  __, 

10.45 

.0005225 

0.96 

i 

TABIE  6.  Search  widths  for  various  intensities. 


66 


r*etf 


o  L 


* fh^^Y"  >/\ly*> 


Fig.  19  Identification  of  symbols  for  previous  tables. 

Interpretation  of  above  tables i 

For  scaling  purposes  use  column  I2g  as  the  surface  intensity  which 
is  known  and  from  which  it  is  desired  to  scale. 

For  example,  if  the  contour  of  a  pool  of  intensity  .05  r/hr  is  known 
and  it  is  desired  to  determine  the  search  width,  enter  column  L^q  and 
find  .05  r/hr.  The  search  width  would  be  the  pool  size  plus  .6122  miles 
when  searching  at  500  feet  altitude  or  plus  1.1 8  miles  (interpolating) 
for  a  search  height  of  1000  feet. 
Conclusion 

The  diameter  of  the  detector  view  can  be  added  to  the  pool  diameter 
to  obtain  search  width  when  scaling  from  surface  radiation  more  than  one 
order  of  magnitude  larger  than  the  minimum  circular  surface  intensity. 
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APPENDIX  IV 
EFFECT  OF  OCEAN  CURRENTS 
A.  General. 

There  are  many  contributing  causes  of  currents,  both  external  and 
internal  to  the  water.  In  the  deep  open  ocean  wind  is  a  prominent  in- 
fluence. Persistent  winds  interact  with  surface  water,  imparting  up  to 
two  percent  of  the  wind  velocity  to  water  movement.  The  direction  of 
the  movement  is  complicated,  with  the  resultant  of  the  "coriolls  force" 
and  layer  friction  balancing  the  actual  wind  imparted  force. 

P.  H.  Kuener  in  MARINE  GEOLOGY  /V,  page  29,  is  quoted  in  part: 

Flow  in  deep  water  far  from  the  coast  driven  by  a  constant  wind 
(i.e.,  drift  currents)  at  a  constant!  speed,  the  motion  at  the  sur- 
face is  directed  4.5°  to  the  right  (in  the  northern  hemisphere)  of 
the  wind.  Below  the  surface  the  direction  gradually  diverges 
farther  and  farther  to  the  right  of  the  wind  and  concomitantly  the 
velocity  decreases.  At  a  depth  called  the  friction  depth,  which 
varies  with  the  latitude  and  force  of  the  wind  and  does  not  exceed 
200  m,  the  current  is  directly  opposite  to  the  surface  current  and 
has  fallen  to  L$   of  the  surface  velocity.  If  bottom  friction  is 
left  out  of  account  the  average  movement  is  at  right  angles  to  the 
direction  of  the  wind.  The  frlctlonal  force  exerted  by  the  wind 
and  the  coriolls  force  at  right  angles  to  the  movement  are  then 
equal  and  exactly  opposed  in  direction.  This  state  is  closely 
attained  in  nature.  Yet  there  is  some  slight  friction,  and  the 
resulting  force  plus  the  "coriolls  force"  is  equal  and  opposite  to 
the  force  exerted  by  the  wind 

When  the  wind  just  begins  to  blow  the  first  movement  of  the  water 
will  be  in  the  same  direction  as  the  wind,  but  as  soon  as  motion 
of  the  water  sets  in  the  rotation  of  the  earth  makes  itself  felt. 
In  shallow  water,  where  the  friction  with  the  bottom  is  consider- 
able, the  current  will  be  less  deflected  by  the  coriolls  force  than 
in  the  ideal  ease.  In  regions  like  the  North  Sea  and  the  Baltic, 
where  the  winds  are  variable,  where  the  water  is  shallow,  and  where 
the  surrounding  coasts  are  at  no  great  distance  wind  driven  currents 
will  generally  follow  the  direction  of  the  motive  force  fairly 
closely. 

Internal  characteristics  and  variants  create  an  isobarlo  gradient 

between  waters.  Salinity,  temperature,  pressure,  level,  etc.  all  have  a 

bearing  on  currents;  however,  considering  all  open  ocean  currents,  none 
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are  stronger  than  the  GULF  STREAM  which  has  a  velocity  of  up  to  3  knots. 

Most  open  water  movements  run  between  0.3  and  0.7  knots,  but  it  is  very 

difficult  to  predict  the  direction  of  flow  at  any  given  location  at  any 

given  time. 

Bowditoh  /kl  states: 

Current  arrows  on  nautical  charts  represent  average  conditions  and 
should  not  be  considered  reliable  predictions  of  the  conditions  to 
be  encountered  at  any  given  time. 

8.  Analysis  of  Currents  in  Sample  Ocean  Areas. 

1.  Set. 

Enclosed  are  data  on  current  direction  and  speed.  The  first 
eight  (8)  pages  (figures  20  -  27)  are  plots  of  the  frequency  of  current 
directions  (set)  for  the  specified  areas  for  the  months  indicated.  Data 
were  collected  by  the  Hydrographic  Offioe  from  ship  reports.  The  number 
of  observations  upon  which  the  graphs  are  based  is  indicated  inside  each 
curve.  This  figure  is  followed  by  a  single  number  representing  the  per- 
cent of  observations  which  reported  no  current  (no  set  or  drift).  The 
area  represented,  broken  into  5°  X  5°  quadrangles,  is  a  ten  degree  strip 
of  latitude  (25°N  to  35°N)  across  the  Pacific  from  150°E  to  170°W  longi- 
tude (about  1-£-  million  square  miles).  Choice  of  the  area  was  Influenced 
by  the  large  number  of  observations  available  over  such  a  large  expanse 
of  open  ocean.  The  months  plotted  were  chosen  arbitrarily  to  provide  a 
three  month  period  to  observe  gradual  changes  (January,  February,  and 
March)  and  other  samples  during  the  year  to  show  overall  variation  (June 
and  September).  The  source  of  the  information  is  the  ATLAS  OF  SURFACE 
CURRENTS  /5/. 

Such  random  behavior  lends  itself  to  a  probablistic  formulation 
rather  than  a  deterministic  approach.  However  the  metamorphosis  of  the 
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frequency  plots  precludes  any  distribution  other  than  UNIFORM.  This 
assumption  is  the  safest  one  to  make  when  generalising  for  large  areas 
and  over  the  whole  year, 

A  bar  chart  of  current  persistency  derived  from  Hydrographic  records 
and  compiled  by  Mrs.  Martha  Olson  of  USNRDL,  San  Francisco,  Calif,  is 
enclosed  as  figure  28.  This  chart  shows  the  frequency  in  large  ocean 
areas  that  the  most  frequent  set  in  smaller  5°  X  5°  quadrangles  occurs 
for  three  categories  of  persistency.  For  example:  In  55%  of  the  5°  X  5° 
sub  areas  of  the  Northwestern  Pacific  the  most  frequent  direction  of  the 
set  occurs  in  less  than  25/6  of  the  observations.  In  4-5/6  of  the  sub  areas 
the  most  frequent  direction  of  the  set  occurred  in  25/6  to  50%  of  the  ob- 
servations. Table  7  lists  the  data  from  which  the  bar  chart  was  derived. 

2.  Drift. 

Figure  29  was  compiled  by  Mrs.  Olson  from  1961  Hydrographic 
Office  information  not  yet  published.  It  depicts  the  frequency  of  drift 
for  various  mean  speeds  of  the  "most  frequent  current. ■  It  shows  that 
in  over  80%   of  the  areas  the  drift  of  the  most  frequent  current  is  be- 
tween 0.4  and  0.8  knots.  This  graph  refers  only  to  the  North  Atlantic. 

Table  8  that  follows  lists  the  mean  speed  and  deviation  for  various 
sectors  at  various  times  of  the  year.  In  general  the  deviation  is  no 
more  than  0.2  knots,  and  the  mean  is  less  than  the  mean  of  the  speed  in 
the  most  frequent  direction  (because  when  the  set  is  in  other  than  the 
most  frequent  direction  there  is  a  tendency  for  the  movement  to  be  less 
defined. 
C •  Conclusions .  , 

1 .  SET  is  an  unpredictable  variable  at  any  given  time.  The  safest 
and  broadest  assumption  that  can  be  made  is  that  the  direction  of  current 
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at  any  given  time  and  at  any  given  place  is  UNIFORMLY  die tributed  . 

2.  The  strip  considered  in  figures  20  through  27  was  sampled  for 
the  month  of  January  to  determine  the  percent  of  area  under  the  curve 
within  the  semicircle  centered  at  the  prevailing  set.  Two  thirds  of  the 

total  area  is  included  under  this  half  of  the  curve. 

i 

3.  No  matter  what  the  direction  of  the  set  is,  the  drift  is 

reasonably  predictable.  The  variance  of  the  drift  is  small  and  in  view 
of  the  time  late  involved  in  the  problem,  the  solution  is  not  sensitive 
to  such  small  variances. 


71 


£Uu*r<juj5l-        «u.      OcQ^ui<^N,.ctK-~o~ivT 


72 


ClU^OW^^       Ok       o  00  MUi  ^/  ^^t  V--0^^ 


Q-tu^o**^^     ou_      Oofl^^^aK-OZ^) 


^**/«tf  oik,  ^.  i- 


u_      o  <£>  ia  **>  <XL  "i>^  *-    -  O  ?  vO 


O^Wl><*  \JLJ  ~Z\-  O  U.       O  Qa  \^tii  ^    ^»  «fc  I  -    ^  ,r   5. 


7* 


A-  VjjQf  u  ui  ^  |^  <$i± 


i 


MM 


3    oO  vo**1  ^  "*  <  *-  -  S* 


07 


"1 


8* ! 


PP1 


^fi^;i 


77 


0-W<*C-vjUJ  2?!"         OU       0<Dnufci/><tl o£^ 


2-  Uj  Q£  o  U.,  Z  |-        OH        Od)tiU  e/>«t^ 


>.  w-i 


■ 


■nzd.co^eA/T  p£«s/sr£tfcy    t*w"1  jKtf^&gw^ 


&L    ojl  Qbserytkfioy 


iS     tr> 


uch  rvo$t   prc%  ^ ?yi,T  direc?/ o v\  occurrea 


i 


TABLE  7 

A  MEASURE  OF  CURRENT  PERSISTENCY  IN  FIVE  SELECTED  AREAS 

PART  1 

Number  of         Number  of  Number  of  Quadrangle-Months 

Quadrangle-    Observations     in  Which  Most  Frequent  Di- 

REGION  Months  rection  is  Observed  the 

Indicated  *  of  Observations 

25*         25*-5Q*  50* 


9,512 

12,595 

184,760 

1,092 

10,236 


A. 

Northwestern 
Pacific  Ocean 

80 

B. 

Northeastern 
Pacific  Ocean 

150 

C. 

North  Atlantic 
Ocean 

253 

D. 

East  Indian 
Ocean 

59 

E. 

Northwestern 
Indian  Ocean 

96 

45 

35 

0 

62 

78 

10 

90 

149 

14 

27 

Jl 

1 

28 

64 

4 

TOTAL 


638 


218,195 
PART  2 


252 


357 


29 


*  of  Quadrangle-Months  in  Whioh  Most  Frequent  Direc- 
REGION                     tion  is  Observed  the  Indicated  *  of  Observations 
2#  25fr-5(fl  50? 

A.  Northwestern 
Pacific  Ocean 

B.  Northeastern 
Pacific  Ocean 

C.  North  Atlantic 
Ocean 

D.  East  Indian 
Ocean 

E.  Northwestern 
Indian  Ocean 


56* 

44* 

0* 

41* 

52* 

7* 

35* 

59* 

6* 

66* 

29* 

5* 

29* 

67* 

4* 

TOTAL  39*         56* 

♦Additional  information  on  the  following  page. 
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TABIE  7  -  Cont«d. 

Quadrangle-Month i  Grouping  of  data  from  all  observation*  over  tha 
years  taken  during  the  indicated  calendar  month  and  in  the  sane  quad- 
rangle of  approximately  5°X5°. 

Source  of  Data:  Measurement  of  lengths  and  widths  of  spokes  of 
current  roses  in  HO  Publications,  Atlas  of  Surface  Currents, 


Courtesy  of  Mrs,  Martha  Olson, 

TJ5NHDL,  San  Francisco,  Calif. 
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APPENDIX  V 
SEARCH  PLAN  FOR  SMALL  FIX  DEVIATION 
Two  other  errors  in  location  are  considered  besides  the  fix  error. 

a.  Navigational  errors  of  an  aircraft  in  finding  the  exact 
location  to  which  dispatched  are  characterised  by  a  normal  distribution. 
In  this  model  it  is  assumed  that  the  error  due  to  navigation  is  normally 
distributed  with  a  standard  deviation  of  five  miles.  The  variances  of 
the  fix  error  and  navigation  error  are  additive  with  a  resulting  standard 
deviation  of  5.6  miles  C||(2.5)2+(5)2). 

b.  Current  flow  is  assumed  to  be  equally  likely  in  direction 
at  a  given  location  and  time.  The  drift  is  dependable  enough  as  pub- 
lished for  the  area  of  interest. 

The  first  pass  of  the  area  should  conform  to  the  pattern  shown  in 

figure  31 •  The  search  is  designed  to  investigate  the  peak  of  the  dis- 

i 
tribution  at  the  time  of  arrival  and  follow  up  with  loops  outside  and 

then  inside.  This  will  maximise  the  probability  of  detection  per  time 

of  search  effort.  Koopman  /2/  discusses  this  plan  in  detail, 

u  =  drift 

v  =  aircraft  velocity 

T  =  time  elapsed  since  initial  detection 

The  time  to  go  from  A  to  B  in  Figure  30  is  r-j-t-^ 

T 
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L* 


L 


7 


Fig.  30  Pattern  for  retiring  search. 

The  time  for  the  object  of  the  search  to  more  from  circle  r-j  to 

circle  r2  =  r2"r1 
u 

To  make  this  equal  to  search  vehicle's  time  to  go  from  A  to  B, 

rg  =  mr-j,  where  m  ■  v+u.  Similarly  r^  =  mr2J  *%  =  mroj  etc.... 

v-u 

To  prevent  overlap  and  expand  or  contract  about  the  loop  re,  /  mr/. 

To  make  this  leg  be  S  miles  outside  the  first  leg  determine  re  such  that 

rg+r^  =  r*r*7-S»  from  which  re,  =  mr/+a  where  a  =yfl 
V      u  ^^ 

The  distance  r^  should  equal  uT.  However  since  in  practice  the  ex- 
panding square  is  used  to  approximate  a  spiral,  a  mean  r-j  is  selected  so 
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that  r-j  a  0.9ut.  To  further  allow  for  the  fact  that  any  change  in  course 
of  the  pool,  at  the  same  speed  will  reduce  the  outward  component,  r-j ,  it 
selected  as  0.8ut. 


L*  =  XBT-i 

l>2  -  mLi+r-j 

h  =mI2 

L/  =  mLa+a 

L5  =  "fy 


(outside) 


Iv  =  mlr4a 

Lg  =  mLy-2a   (inside) 

1^  =  mLg 

L  q  =  ml«— 2a 

L11  =  mL,0 

L|2  =  mLj.j+3a  (outside) 

etc. 


A  sample  search  is  worked  out  for  the  following  situation: 
Height  of  aircraft    =  500  feet 
Spacing  (S)         =  9.61  miles 
Drift  (u)  =  0.5  knots 

Speed  of  aircraft  (r)  =  180  knots 
Time  late  (T)        =24-  hours 
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r-,  =  .8(.5)(24)  =  9,6  miles 

a    =  yj  =  I80(9.frl )  =  9.65 
T-a  179.5 

■    =  I±B  =  160.5  =  1 .0056 
v-u      179.5 

1*1  =9.65  Total  mileage  for  this  pattern  460  mi. 

Ig  =  19.3  Time  for  pattern  -  2  hrs  35  min. 

1*3  =19-4  If  this  plane  only  has  3  hrs  on  station 

L^  =  29.2  available  it  should  spend  the  last  25 

L5  =29.4  minutes  sweeping  the  lower  left  corner. 

I«£  =39.3  Another  aircraft  will  have  to  make  the 

1*7  =  39.5  parallel  sweep. second  pass  which  re- 

Lg  =  20.4  quires  1  hr  40  min  (300  miles). 

1^  ■  20.5 

L11  =  n331 
1,2  =  30.3 
L,3  =  30.5 
I14  =  59.7 
1^5=60 
L,6  (end)  =  60 

The  plan  leaves  a  holiday  in  the  lower  left  corner  of  the  pattern. 
In  order  to  take  advantage  of  the  best  information  the  upper  right  corner 
should  be  oriented  in  the  direction  of  the  prevailing  current  so  that  the 
starting  leg  is  45°  to  the  left  of  the  most  frequent  current. 

Note  that  the  oenter  may  be  adequately  oovered  before  outside  areas 
are  searched  -  depending  on  the  balance  of  the  time  late  and  standard 
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deviation  of  location.  In  ease  this  happens  the  pattern  should  be  con- 
tinued with  outside  expansion  only.  In  the  example  cited  above  the  cen- 
ter search  is  just  a  narrow  turn  of  270°  to  right. 

After  this  search  is  completed  proceed  with  a  parallel  sweep  as  ex- 
plained in  paragraph  11b  of  the  main  body  of  this  st\Jdy. 
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APPENDIX  VI 
SEARCH  PLAN  FOR  LARGE  FIX  DEVIATION 
For  the  case  where  the  fix  Is  not  any  more  certain  than  within  a 
standard  deviation  of  25  miles  or  more,  the  errors  in  navigation  and 
current  movement  can  be  ignored.  A  set  of  expanding  squares  -  squares 
of  uniform  coverage  -  are  superimposed  in  the  square  pattern.  Koopman 
/2/  covers  the  subject  in  detail. 

The  altitude  of  the  aircraft  should  be  as  low  as  operationally 
acceptable  and  for  an  example  500  feet  will  be  used. 

Figure  32  is  a  drawing  of  the  pattern  for  the  following  problems 
Aircraft  speed  =180  kts. 
Endurance  on  station  >  54-0  miles 
Spacing  =  11  .0  miles 
Average  search  width  (W=E)  =  8.44.  miles 
Coverage  factor  0  =  WL  =  4558  square  miles 


Maximum  search  radius  =  a  ■  (43.7  miles  =  M  '         !  )  i 


■  t&> 


Half  length  of  Ejj»  square  (S^)  = 

Number  of  passes  n  =  a  S  =  2 

2pW 

51  =  19.5 

52  =  33.8 
Length  of  a  side  of  K^  square  =  2Sfc 

1^   =  39 
Lg  =  67.6 


r-tr  v   i  <-"i- 


V; 


-;]  *■ 
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Number  of  circuits  vithin  K*h  major  square  (N^)  =  Lk 

2S 

N-,  =  2 

N2  =  3 

The  second  uniform  square  is  tilted  U5°  to  the  first.  This  pattern 

is  plotted  in  figure  32. 

The  average  probability  of  detection  ist 

P  =  1-(l+nw]eocp/^aWl=  .&U 
o     *^  o 

If  the  search  aircraft  has  more  endurance  than  considered  above  (3 

hours  at  180  kts  -  assuming  the  search  area  up  to  1000  miles  from  base) 

the  plan  should  be  appropriately  recomputed. 

For  instance  if  search  area  were  only  500  miles  from  base 

Station  Endurance   1500 

Coverage  Factor  0  =  WL  =  12660 

Maximum  search  area  radius  a  =  55.4- 

Number  of  passes  n  =  3 

51  =  19.5  Li  =  39  R,  =  2 

52  =  33.8  Lj  =  67.6  N2  =  3 

53  =  -43.6  Lj  =  87.2  N3  =  U 
The  average  probability  in  this  case  would  be: 

P  =  1-(l+2wWpp2jfj  =  0.67 
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